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Abstract 


An analytical/experimental investigation was performed to study the effect of material nonlin- 
earities on the response of composite tubes subjected to combined axial and torsional load- 
ing. The effect of residual stresses on subsequent mechanical response was included in the 
investigation. Experiments were performed on P75/934 graphite/epoxy tubes with a stacking 
sequence of [15/0/ ± 10/0/ —15], using pure torsion and combined axiai/torsional loading. In 
the presence of residual stresses, the analytical model predicted a reduction in the initial 
shear modulus. Experimentally, coupling between axial loading and shear strain was ob- 
served in laminated tubes under combined loading. This phenomenon was predicted by the 
nonlinear analytical model. The experimentally observed linear limit of the global shear re- 
sponse was found to correspond to the analytically predicted first ply failure. Further, the 
failure of the tubes was found to be path dependent above a critical load level. 
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1.0 Introduction 


The use of composite materials in tubular structures has become popular in recent years. 
Tubes made from graphite/epoxy have been used in the INTELSAT 6 Communication Satellite, 
and they are also under consideration by NASA for use in the truss structure of the space 
station 1 . In addition to spacecraft applications, composite tubes are also being considered for 
automobile drive shafts and for robotic arms. For applications such as these where high 
stiffness and dimensional stability are critical, a thorough understanding of the response of the 
composite tubes subjected to combined thermomechanical loading is imperative. 

Polymer-based materials, such as the graphite/epoxy used in these tubes, frequently 
exhibit nonlinear constitutive behavior of the reversible or irreversible kind, depending on the 
direction of loading. This nonlinearity can have a significant effect on the distribution of 
stresses and strains within a tube. For this reason, an analytical tool employing a sufficiently 
general nonlinear constitutive theory is desired for predicting the response of composite tubes 
subjected to combined mechanical loading. This kind of tool will allow for a more accurate 
prediction of the state of stress, strain, and deformation in a composite tube. 

Load path dependence is a phenomenon which occurs in nonlinear materials with 
dissipative response characteristics. This is defined by the state of stress and strain at a 
given state of loading being different for each different path taken to reach that state of load- 
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ing. An investigation of this phenomenon at the gioba! level will be possible with an appro- 
priate nonlinear analytical model at the ply level. 


1.1 Literature Review 


Much work has been done in recent years to develop analytical models that predict the 
response of anisotropic cylinders. Many different approaches such as shell theory, finite el- 
ement analysis, and exact and variational approaches to elasticity theory have been utilized. 
Different variations of shell theory have been used to investigate the response of composite 
cylinders to torsion, tension {or compression), bending or a combination of these 2 ” 8 . These 
analyses have yielded generally good predictions for thin wailed cylinders (r 0 /t >10) where 
through-the-thickness strains can be neglected. 

Rizzo and Vicario 9 used finite element analysis to predict the response of laminated tubes 
consisting of generally anisotropic layers with all coupling terms included. Unlike the shell 
theories described previously, this analysis included the effects of gripping the tube. The au- 
thors investigated the effects of varying the thickness-to-diameter and length-to-diameter ra- 
tios for different orientations of unidirectional tubes. In addition, they investigated the effect 
of different gripping arrangements on the subsequent mechanical response. 

More recently, elasticity solutions for the composite tube problem have been developed. 
Cohen and Hyer 10 developed three axisymmetric elasticity solutions to determine residuai 
stresses in cross-ply tubes: plane stress, plane strain, and generalized plane strain. They 
found that the generalized plane strain solution provided the best predictions of the three. 
The results from the elasticity solutions were compared with results from classical lamination 
theory, and it was found that there are limitations to approximating the response in a tube 
with a flat plate solution. 


Introduction 


2 



Hyer, Cooper, and Cohen 11 subsequently developed a solution to predict the response 
due to uniform temperature change based on the principle of complimentary virtual work. The 
results from this solution were compared with the elasticity solution developed previously 10 
and agreed quite well. This work was restricted to cross-ply tubes. 

The elasticity solution for cross-ply tubes was extended by Hyer and Cooper 12 to investi- 
gate the problem of a cross-ply composite tube subjected to a circumferential temperature 
gradient. Finally, the axisymmetric elasticity solution was extended by Rousseau, Hyer, and 
Tompkins 13 to predict the response of arbitrary laminated tubes subjected to mechanical loads 
as well as uniform temperature changes. 

All of the analytical tools described above have employed a linear constitutive theory. 
Orgill and Wilson 14 have developed an analytical model to approximate the response of tubes 
using a nonlinear analysis. A strength of materials approach was used for analyzing tubes 
consisting of a single orthotropic layer. Both geometric and material nonlinearities were 
considered. The geometric nonlinearities were handled using an incremental procedure. At 
each load increment, the finite changes in the tube geometry { wall thickness, inside diameter, 
outside diameter, etc. ) were computed and adjustments made to the cylinder dimensions. 
The material nonlinearities that were considered arose from finite changes in the orientation 
of the fibers in the layers due to deformations in the tube. The nonlinear constitutive behavior 
was determined from experimental methods and approximated by 4- Be^ + Cs^ with 

similar expressions for a 22 and This type of approximation assumes that there is no 
interaction between the stress components and that the stress-strain response is elastic. 

Loading path dependence on the response of composite materials has not been so thor- 
oughly investigated. Loading path dependence as it relates to the stress state at failure of 
composite tubes has been investigated by Choo 15 . This research was limited to hoop wound 
( 90° ) tubes. An extensive experimental program was carried out in which tubes were tested 
in combinations of compression and torsion or compression and internal pressure. It was 
found that in the fiber fracture mode, the state of failure stress was load path independent. 
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In the matrix fracture mode the state of failure stress was found to be load path Independent 


in some cases. 


1.2 Objectives 


The objective of the present investigation is to characterize the response of composite 
tubes under combined mechanical loading. Towards this end, and analytical/experimental 
program will be undertaken. The analytical program entails the development of a model ca- 
pable of predicting the response of composite tubes subjected to combined thermomechanical 
loading. The analytical model will be developed In a general fashion to include different types 
of mechanical loading such as axial force, torque and internal and external pressures. Ther- 
mal effects will also be included to account for residual stresses induced by the fabrication 
process. An elasticity approach based on the solution obtained by Rousseau, et al., 13 will be 
used. The present model modifies this solution with the incorporation of a nonlinear 
constitutive theory, in addition, an efficient algorithm based on the structural matrix approach 
was incorporated for solving the resulting equations. This algorithm is particularly useful for 
laminated tubes with a large number of layers. 

The predictions of the model will be correlated with the results of an experimental pro- 
gram. The experimental program will specifically investigate the response of composite tubes 
subjected to combined axial and torsional loading. The tubes used in this investigation are 
two inch nominal diameter tubes with a wall thickness of approximately 0.060 inches. They 
are fabricated using P75/934 graphite epoxy tape with a stacking sequence of 
[15/0/ ± 10/0 — 15],. This configuration is one under consideration by NASA for use in the 
truss structure of the Space Station. 
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2.0 Analytical Model 


2.1 Introduction 


Rousseau, et al., 13 developed an exact planar elasticity solution using a linear constitutive 
theory to study thermally induced stresses and deformations in composite tubes. In the 
present investigation, the effects of material nonlinearities will be studied. Towards this end, 
a general nonlinear constitutive model will be incorporated into the solution obtained by 
Rousseau, et al., to permit the usage of any nonlinear constitutive theory that is desired. Once 
the nonlinear solution is complete, a suitable nonlinear model will be incorporated. 

Due to the nonlinear effects, a modification to the numerical solution technique is desired 
to improve the computational efficiency of the analytical model. For this purpose, the 
Local/Global Stiffness Matrix formulation will be incorporated. Although thermal effects and 
the effects of internal and external pressures will not be included in the experimental portion 
of this work, the analytical formulation will Include the presence of these types of loading for 
more generality. The material properties will be assumed to be temperature independent. 

Figure 1 shows the coordinate system and nomenclature used for a laminated tube. The 
inside and outside radii of the tube are denoted by r„ and r N respectively, where N is the total 
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number of layers In the tube. The radii of the interfaces between layers are denoted as r,. r 2 , 
etc. The governing equations to be derived will be developed in cylindrical coordinates; x 
represents the axial coordinate, 6 represents the circumferential coordinate, and r represents 
the radial coordinate. The displacements in the x, 6, and r directions will be represented by 
u, v, and w respectively. The ply orientation angle is <f>, and the principal material coordinates 
are oriented such that the 1-axis corresponds to the fiber direction, the 2-axis is normal to the 
fiber but in the lamina surface, and the 3-direction is normal to the lamina surface. The 
internal and external pressures are denoted by P, and P 0 respectively, and the applied axial 
and torsional loads are denoted as F x and T x respectively. 


2.2 Formulation 


2.2.1 Field Equations 


In the equations that follow, the contracted notation for the stresses and strains in the 
principal material coordinates will be used, l.e. 


- 


- 


• 


" 

*1 
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®i 


e ii 
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713 
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For the purposes of this derivation, the field equations will be considered for a single 
layer and then combined to form a solution for a complete laminated tube. The equations to 
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be used are the Equilibrium Equations, Compatibility Equations, Strain-Displacement 
Equations, and Constitutive Equations. The Equilibrium Equations in cylindrical coordinates 
are 


do r 

dr 


+ -p (<r r - o d ) -f 


1 dj xr 

r de + ax 


+ 


R r = 0 


[2.1a] 


dr 9r 1 gffl ^ , 2 , D 

ar + r aa + ax + r T *r + - 0 


[2.1b] 


d T xr 1 d<x x 

ar + r aa + ax 


+ 


k xr 


+ R X = 0 


[2.1c] 


where R x , R fl , and R r are the respective body forces. The Strain-Displacement Equations are 


y&r 
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£x ax 
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£fl = "r'('i' + w ) 

[2.2b] 

aw 
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The general form of the Compatibility Equations in cylindrical coordinates are 


*r , A x _ g2 Vxr 
ax 2 aV dxdr 
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[2.3a] 
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[2.3b] 
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The constitutive equations in the principal material coordinates are: 
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or the inverse 
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[2.4b] 


where C (J are the initial stiffnesses of the lamina, are the initial compliances, sj are the 
strains induced by a change in temperature, and the sf* L are the nonlinear strains. This 
equation can be transformed to the global cylindrical coordinates of the tube. The result is 
as follows: 
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where C, } are the transformed initial stiffnesses of the lamina. This equation can be inverted 
to give strains in terms of stresses. 
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where Ay are the transformed initial compliances of the lamina. For the moment, the nonlinear 
forms of the constitutive equations will be left in a general form with no specific nonlinear 
model used. This will be incorporated later in this chapter. 


2.2.2 Boundary and Continuity Conditions 

The field equations will be used to determine a general form of the governing equations 
for the displacements in the tube. To obtain a specific solution, the boundary conditions will 
be applied to the governing equations. On the outside surface of the tube (r = r N ), the tractions 
will be specified as follows: 


Y(x. 0. I'm) == - p o 

[2.6a] 

*er(*. 0 * r N)~O 

[2.6b] 

TxrCx, 0 , r N ) = 0 

[2.6c] 


Similarly, the tractions on the inside surface (r = r 0 ) will be specified as 


a r (x, 0 , r 0 ) = - P, 

[2.7a] 

T 0r (x. e > r o) = 0 

[2.7b] 

r xr (x, 0 , r 0 ) — 0 

[2.7c] 


For an arbitrary N layered tube, continuity of tractions and displacements at the interfaces 
of the layers must be fulfilled. The continuity of tractions can be expressed as 

<4 k) (x, 0. r k ) = <r< k+1) (x. e, r k ) k = 1 ,2 N - 1 [2.8a] 

t^x, 6, r^) = T^ r +1) (x, 6 , r k ) k = 1,2 N-1 [2.8b] 
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k— 1,2 N — 1 


[2.8c] 


r^(x, 0 , r k ) — T xr ^ ) (x, 0 , r k ) 
and the continuity of displacements are given by 


u (k) (x, 0 , r k ) = 

= u (k+1) (x, 9 , r k ) 

k— 1,2,.. 

.,N-1 

[2.9a] 

v< k >(x, 6 , r k ) = 

= v (k+1) (x, 9 , r k ) 

k — 1,2,.. 

.,N — 1 

[2.9b] 

w (k) (x, 0 , r k ) = 

= w* k+1 *(x, 0 , r k ) 

k= 1,2,. 

...N — 1 

[2.9c] 


where k represents the k th layer in the tube and ranges from 1 to N-1. Finally, the total axial 
force, F„, will be introduced to the tube as 

f2rr rr N 

F x = I a r (x, 0 , r)rdrd# [2.10] 

**0 Jr„ 

and the applied torque, T x is 


n r N 

*x <?( x . 8 ’ r)r*drd0. [2.11] 

0 


These two equations will be referred to as the "integral force and torque equations." 


2.2.3 Assumptions 

The applied internal and external pressures are Independent of 0 and the temperature 
will be assumed to be spatially uniform. This leads to an axisymmetrlc problem, l.e. the 
stresses, strains, and displacements are independent of 9 . For the portion of the tube that is 
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sufficiently far from the point of application of the axial force, F x , and the torque, T x , the 
stresses and strains will also be independent of the axial coordinate. These assumptions can 
be expressed as 



[2.12a] 



[2.12b] 


where 


X = stress, strain, or displacement 
Y = stress, or strain 


2.3 Solution 


2.3.1 Solution for Governing Equation 

The governing equations will be derived for an individual layer. Application of boundary 
and continuity conditions will yield a system of equations for the radial displacements at the 
interfaces of the individual layers. Determination of the interfacial radial displacements in 
terms of the applied boundary conditions is sufficient for determination of the response of the 
tube in terms of global and local parameters. 

Incorporating the assumptions of axisymmetry and the uniformity of the stresses and 
strains in the axial direction (which ensures that the radial displacement is a function of r 
only), the displacements take the form 
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u = u(x,r) 


[2.13a] 


v = v(x,r) 


w = w(r). 


Using these assumptions, the field equations can be simplified as follows: 


Strain-Displacement Equations: 
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dx 
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dw 
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dr 



Compatibility Equations: 



1 *x 
r dr 


= 0 


[2.13b] 

[2.13c] 

[2.14a] 

[2.14b] 

[2.14c] 

[2.14d] 

[2.14e] 

[2.14f] 

[2.15a] 

[2.15b] 
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£[t -fM-o 


[2.15c] 


where the three remaining compatibility equations are satisfied identically. 
Equilibrium Equations: 


~ + “r* ( ff r *“ °e) — 0 

[2.16a] 

d *9r ,2 n 

ar + r T *r~° 

[2.16b] 

dr * r 1 n 

ar + 1* Txr ~° 

[2.16c] 


where the body forces have been neglected. 


These twelve equations ( [2.14], [2.15], and [2.16] ) will now be used to derive the general 
expressions for the displacements of a given lamina. Equation [2.15a] t in conjunction with 
Equation [2.12b], can be integrated to solve for the axial strain, e* directly. 

®x( r ) = Ar + B 

where A and B are constants. Equation [2.15b], however, requires that A=0. Renaming B 
as e 0 , we find that the axial strain within the lamina is a constant. 

e x ( r ) = e o [2.17] 

Using equations [2.14a] and [2.17], the axial displacement u{x,r) can be solved for. 

du 

e * = 17 = 6 ° 

u(x,r) = e 0 x + f(r) [2.18] 
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Equation [2,15c] can be integrated to obtain 

Vxtf(r) = Cr + -B- [2.19] 

as the in-plane lamina shear strain. From equations [2.14f] and [2.19], the tangential dis- 
placement v(x,r) can be obtained. 


v - Jv_ _ rr , SL 

^ ~ Cr + r 


v(x,r) = ^Cr + “-)x + g(r) 


[ 2 . 20 ] 


Equations [2.16b] and [2.16c] can be solved for and t* respectively. This yields 


l xr ““ x 


[ 2 . 21 ] 


T 0r “ 


[ 2 . 22 ] 


Inserting equation [2.18] into [2.14e], we find that 


Vxr 


m 

dr 


From the constitutive relations 


NL 


Vxr — A 45 T 0r + A 55 T xr + Vxr 


substituting equations [2.21] and [2.22] and integrating, we find 


f( r ) — 


A45F 


-F A55E In r + F 1 + 




[2.23] 
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where F t is the rigid body translation. 

Equation [2.20] can be inserted into equation [2.14d] to obtain 

<*9(0 g(r) 2D v 

y 0r~~ dr r r 2 x 


But, since the strains are functions of r only, D = 0, therefore, 

dg(r) g(r) 

v#r “ dr r 


[2.24] 


Using the constitutive relations 

y 0 r = A44T 0r + A45T xr + ^ r L [2.25] 


Equating [2.24] and [2.25], and solving for g(r), we obtain 


fl(r)»- 


A 44 F “ 

2r 


A45E ■+• G<j r + r 


NL 

yor 


dr 


[2.26] 


where G 1 is the rigid body rotation about the x-axis. Substituting these expressions into the 
expressions for the displacements, u and v, the resulting expressions for the displacements 
are 


u(x,r) = e 0 x 


A45F - 


-f A^E In r + 


1 


y!J. L dr+Fi 


[2.27] 


, S A 44 F 

v(x.r) = y Q rx — 


— A45E + r 


J+* 


+ G-jr 


[2.28] 
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where the constant C has been renamed y 0 . The physical meaning of y 0 is the angle of rotation 
of the tube per unit length. Transforming y" L and y$-tothe principle material coordinates and 
neglecting rigid body motion gives 


A 45 F 


h 


u(x,r) = £q x - + A 55 E In r + | ( -ne 4 L + me^dr 


A 44 F I (ms 4 u + nsg L ) 

v(x,r) = y 0 xr - - A^E + r | -- --- — — dr 




[2.29] 


[2.30] 


where 


m = cos <f> 


n = sin <j> 

Consider the remaining equilibrium equation, [2.16a], and the constitutive equations for 
a r and <7*. namely 


C 21 ( e x e x ®x) + C 22( s 0 “ s 0 L — ®x) + C 23(®r ~ e ? L “ s r ) + C 26 (?x0 ~ Yx0 ~~ yI$) 

<r r — C 31 (fi x - e x - e x ) -f - eg 1 - ej) + G^e,. — e^ L — &J) + C 36 (y x0 — yj^ — y x0 ) 

Replacing the total strains, e x , & g , & r , and y with their respective strain displacement ex- 
pressions, and transforming the thermal and nonlinear strains into the principal material co- 
ordinates we obtain 


a e - c i 2 


~ au +c 22 ^+c 


dw , dv 

23 ~rr + <-26 


ax ■ r ■ dr ■ -w dK 
[n 2 C|i(e| NL + e^) + m 2 C i2 (e| NL + ef) + 2mnC 66 Eee ] 


[2.31] 
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[2.32] 


~ 3u . ~ W,~ dw , ~ dy_ r / NL . Tv 

a r - C 13"aJ7 + ^3“ + C 33 fa + ^36 C i3V 6 i + 8 i ) 


where the repeated subscript, i, denotes a summation over the range i~ 1,2,3. At this point, 
it must be noted that two separate cases must be considered. These are the transversely 
isotropic lamina, where C 12 = C 13 , C K = C 33 , and 044 = 2 ( 022 - 023 ), and the monociinic lamina. 
When the constitutive relations are substituted into the equilibrium equation we find, for 
transversely isotropic layers: 


2 d 2 w . „ dw 
^~d^~ +r lF 


C|3 n 2 C I1 

c 22 c 22 


C 22 / C 22 j 


[2.33a] 


and for monoclinic layers 


^ d 2 w 


w = K,c 0 r+ Kj Vo r 2 + r 2 - 3 *- 


(n 2 C (1 4 - m 2 C j2 — 4- S| ) 4- 2mnC 66 e 6 


[2.33b] 


where 


'12 “ H3 


0 2 6 — 2030 


Because the distribution of nonlinear strains through the thickness of the lamina is not known 
a priori, we will assume that they can be approximated by a fifth order polynomial, i.e. 


£| NL = 3| 4* bjf 4* Cjr 2 + d j r 4* ejr 4* fjr 
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For the moment we will assume that the constants a,, b ( , c,. d,, e,, and f, are known constants. 
The reason for this assumption will become apparent later. With this assumption we can solve 
the differential equations for the radial displacement In terms of these constants 
a it b t , c„ d„ e,, and f, . For transversely isotropic layers 

r 2 + r - w ~ C^r + Q 2 r 2 + Q 3 r 3 + Q 4 r 4 + Q 5 r 5 4 - Q 6 r 6 [ 2 . 34 a] 

dr ar 


where 


Qi = (a, + ej)~ 2mn 4?^ a 6 

*-22 Cc 


'22 


Q 2 = _2mn-&S-b 9 + (JS- + ^.\ 

C 22 V C 22 C 22 ) 


2Cj3 + F, 


'66 


Q 3 = ( L ]C| - 2mn c 6 

\ ^*22 / '-'22 


Q 4 4i%±AV-2mn#ld 6 


'22 


'22 


^ ( 4C (3 + Fj \ C 66 

Q 5 = I — = ]e, - 2mn — — d e 


'22 


'22 


Q 6 =fi%±n _ 2mn -&! d6 


'22 


'22 


Fj — Cj 3 — n 2 C| 1 — m 2 C i2 


For the monoclinic layers, the differential equation becomes 


0 7s 

2 d w , „ dw c 


dr" 


■+ r 


22 


dr c 


w 


33 




[2.34b] 
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where 


T is , F l( e i + a l) o mn C 66 

T i - K i e o + — = 2mn-~ — a 6 

'-'33 '-'33 


T 2 = K 2 y 0 + - (C| L + ^ b| - 2mn b e 

C33 C33 


T (2C|3 4- F|) Cqq 

To — = Ci 2mn — — c B 

C33 1 C33 6 


(3C ]3 4* Ft) Cgg 

T !i d 2mn-=^d 6 

C 33 C 33 6 


(4C,3 + F|) C 66 

t 5 = = e j — 2mn-= e 6 

'-'33 '-'33 


(5C | 3 4- F|) C 66 

j e - 2mn -= — 

C 33 C. 


33 


These equations can now be easily solved. The result for transversely isotropic layers is 


w(r) - A 1 r 4- — 4- Pi r In r 4- P 2 r 2 4- P 3 r 3 4- P 4 r 4 4- P 5 r 5 4- P 6 r 6 


[2.35a] 


where A 1 and A 2 are arbitrary constants and 


P, = -4— [F,(a, + ej r ) - 2mnC e6 a 6 ] 


2C 


22 


P 2 _ __1 [(F| + D j3 ^b| — 2mnCggbg] 

3C 22 


P 3 ~ C( F i + 2Cj 3 )C| — 2mnCggCg] 
8C22 
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[(F, + 3C i3 )d| - 2mnC ee d 6 ] 


P 4 = — 4r- 
15C 22 

P 5 — [(F f 4C i3 )ej — 2mnC 66 e 6 ] 

P 6 = ^ F| + 5C '3) f i “ 2mnC 66 f 6 ] 

OOL»22 

For monoclinic layers, the result is 

w (r) = A,r x + A 2 r~ A + G^ + G 2 i* + G 3 r 3 + G 4 r 4 + G,/ + G 6 r 6 


where 



G _ ~ C i3) £ o + + Qj) ~ 2mnC 66 a 6 


G 2 


(C 26 — 2C 36 )y 0 4- (C i3 + F|)b| — 2mnC 66 b 6 
4C 33 — C 22 


Go = 


(Fj 4- 2Cj 3 )Cj — 2mnC 66 c 6 


90 


33 


'22 


(F, + 3C i3 )dj - 2mnC 66 d 6 
1 6033 ^22 


(F, + 4C l3 )d| - 2mnC 66 e 6 
9C 33 — C 22 

(F| + 5C| 3 )f| - 2mnC 66 f 6 
^ 60^ — C 22 


[2.35b] 
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Equations [2.29], [2.30], [2.35a], and [2.35b] are the general solutions for the displacements of 
cylindrical lamina. The solution to a multilayered cylinder can now be found using the solution 
for the single lamina. 


2.3.2 Application of Boundary Conditions 

To fully solve a general N layered tube, the unknown constants in the solutions for the 
displacements must be solved. This gives 6 unknowns for each layer, E, F, e 0 , y 0 A v and A 2 . 
Therefore there are a total of 6N unknowns. These must be determined from the boundary 
and continuity conditions. Restating these conditions: 


— Pi 

4 n Vn)=-p 0 

[2.36a, b] 

£ 

T 

o 

4 N r VN)“0 

[2.37a, b] 

4'rVo) = 0 

o 

II 

"z 

[2.38a, b] 


k = 1,2 N - 1 

[2.39] 

4-Vk) = 4 k r +1> ( r k) 

k— 1,2 N — 1 

[2.40] 

ii 

N 

1. 

/-s . 

ex 

t* 

k= 1,2 N — 1 

[2.41] 

u W (r k ) = u <k+1) (r k ) 

k= 1,2 N-1 

[2.42] 

+ 

II 

> 

k= 1,2 N-1 

[2.43] 

w (k *(r k ) = w <k+1) (r k ) 

k= 1,2 N-1 

[2.44] 


Applying the shear stress-free boundary conditions, [2.37a, b] and [2.38a, b], we find that 
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eW-eM-o 
f (1) = f (N) = o 

from Equations [2.21] and [2.22]. Then from the continuity of shear tractions at the interfaces, 
[2.40] and [2.41], we find 

E (k) = F< k) = 0 k = 2,3 N - 1 

Therefore the T 0r and r,, shear stresses, and, hence, <y 4 and <r 5l are zero through the thickness 
of the tube. The equations for the u and v displacements for the k th layer reduce to 


% 

u W (x,r) = 4 k >x + ( -n^r + 

% 

[2.45] 

f , „(k) NLs 

v (k) (x,r) =s yjj^rx 4- r — ^-dr 

* 

[2.46] 

In a later section, it will shown that ef- and eg 11 - are zero throughout the tube due to the fact that 
<r 4 and a 5 are zero. Therefore 

u (k) (x,r) = 4 k >x 

[2.47] 

v (k) (x,r) = y^rx 

[2.48] 

Applying the continuity of displacements at the interfaces of the layers to the u and v compo- 
nents 
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Therefore the values of e 0 and y 0 are Identical in every layer. This means that the axial strain 
and the rotation per unit length are constant through the thickness of the tube, so 


s (k) = e 

6 0 ~ E 0 


(k) 

yo-yo 


Then 

u^(x,r) - e c x [2.49] 

v^(x,r) = y 0 rx [2.50] 

The remaining unknown constants are A\ k \ A£>, e 01 and y 0 , or 2N + 2 constants. The remaining 
boundary conditions to be satisfied are equations [2.36a,b], [2.39], and [2.44]. These give 2N 
equations so two additional conditions are required. These are the integral force and torque 
equations 


n r N 

a x (r)rdrd0 

■ 0 


( *2n 


Tx = 




r x6»( r )r 2 drd0 


o •Y, 


or, written in terms of the stresses in the individual layers 
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4 k) (r)rdr 


[2.51] 


F 


x 



N 

v* r* 

T x = 2^> T<g(r)i*dr [2.52] 

bt J ^ 


There are now a sufficient number of boundary conditions for a solution to be found. 


2.3.3 Local / Global Stiffness Matrix Formulation 


Rather than solving the 2N + 2 equations by substituting the expression for the stresses 
and displacements into equations [2.36a,b], [2.39], [2.44], [2.51], and [2.52] and solving the 
resulting equations for Af>, A£ k >, s 0 , and y 0 , we will manipulate these equations with the objec- 
tive of reducing the actual number of equations to be solved 18 . For simplicity define 

w j[ = w (k) (r k ) w k = w< k) (r k _i) 

-r-Fto 4 k) ~ 


Writing the expression for w£ and w* for transversely isotropic layers 


w k - Afr k + Tf - + P?r k In r k + + P$r£ + P^ + pgrg 


[2.53a] 


A k 


w k = + 7^7 + P W-i ln r k-i + P&k-i + P&k-i + + P^-i + P6 r k-i [2.54a] 
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and for monoclinic layers 


wit = + ~j~+ G?r k + + G^rg + Gjrj + G^r k r 5 + G£r£ 

4 k 


w k = A l*k-i + j 2 4* <3^r k _ 1 + G2 r k-1 + G 3 r k-1 + G 4 r k-1 + G 5 r k-1 + c 
•kii 


Solving these equations for A$ and A$, we find 


R k Rk 


A k __ r k-i r k 
A i ■“ 


detb 


k r k R k ” R k r k-1 


a; = 


det L 


where 


r k r k < 
det k = 7- 1 

K r k-i r k 


r? = - pW m r k - py k - py k - py k - pr k - pk 


k_2 nk_3 n k 4 n k 5 n k 6 


R k ” w k “ p 1 r k-1 ,n r k— 1 ~ ^2^-1 — R3 r k~1 ” ^4 r k— 1 “ R5 r k-1 — *Vk-1 


for transversely isotropic layers and 


k Rk r k-1 k Rk r k 


AJ- 


detL 


A 2 = 


^Rk-R^ii 


deti, 


where 
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; 6 r k-i [2.54b] 



R 4“ ...4’ pk p k 2 pk 3 pk 4 pk 5 pk 6 

k = w k - G 1 r k - G 2 r k - G 3 r k - G 4 r k - G s r k - G 6 r k 

R k = w k - Gi** k _-i - G 2 r k -1 ~ G 3 r k-1 “ G 4 r k-1 “ G 5 r k-1 “ G 6 r k -1 

for monoclinic layers. Expressing the normal radial stress In terms of the displacements using 
equations [2.14], 12.32], [2.49], and [2.50], we obtain 

<4(r) = cfjSo + C-j 3 + C 33 + C§ 6 y 0 r - (H? 3 S? L + H^ L + H^ L + C^ 1 ) [2.55] 

where 


H 1( == C 1t m + C 2 jn "I- 2mnC0| 

H 2i = C 1t n 2 + C 2i m 2 - 2mnC 6i 
H 6i = mn(C 2i - C 1( ) + (m 2 - n 2 )c e , 


If the expressions for Af and A£ are substituted into the expressions for the radial displace- 
ment, w, which in turn is substituted into [2.55], we obtain an expression for the radial stress 
in the k th layer in terms the radial displacements at the inner and outer radii of the layer. From 
this expression, we can easily write expressions for the interfacial tractions for each layer in 
terms of w*, wf, e°, and y 0 . In matrix form, these are 



K 


11 


*21 


u-k i>k i/k 

*12 *13 *14 

i/-k ixk w-k 

*22 *23 *24 



s O 

Vo 


f nl+ f t“ 

I rk+ 

h NL+ h T 


[2.56] 
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where Fn L and Fft. contain alt of the terms that contain the nonlinear constants 
a„ b,, c„ d„ e„ and f, and Ff and Ff contain all of the terms that contain cj. K|J is called the local 
stiffness matrix for a given layer, k. The expressions for these terms are given in Appendix 
A. 


2.3.4 Assembly of Global Stiffness Matrix 


The local stiffness matrices can be assembled into a global stiffness matrix by writing the 
traction equations at the inside surface, the interfaces between each layer, and the outside 
surface. For the inside surface 

«- = Pj = K^Wq 4* K^ 2 w 1 4- Ki 3 e 0 4- K} 4 y 0 — fJjl — f| 

For the k th interface 


<*r + = 4- K 22 w k + K 23 c 0 4- K 24 y 0 - F^t ~ Ft + 

(k+1)— x x x p(k+1)— p(k+1)— 

-c r —Kn W k + K 12 W k+1 + K 13 e 0 +K 14 y 0 -r N l h T 

For the outside surface 

<*4 + = - P 0 = K 21 W N-1 4- K 22 W N 4- K 23 e 0 + K 24 y 0 - fJJJ - Fj + 

The three resulting forms of the traction equations are: 


for the inside surface 


K^w 0 4* K^W! 4- K] 3 e 0 4- K^ 4 y 0 = P ( 4- Fj 4- fJjl [2.57] 


for the k* h interface 
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[2.58] 


^21 w k-1 + ( K 22 + + KifJ ^ k+1 + (K23 + K13 ^ Vo 

+ (K24 + K ?JVo - (Ft + + Fp +1 h + (Fnl + F ( NL 1>_ ) 

and for the outside surface 

^21 W N-1 + ^22 W N + ^23 e o + ^24^0 = — P 0 + F j + + F NL [2-59] 

where 

W k = w k = w k+1 
a^~4 k+1 >- = 0 

In addition to these equations the integral force and torque equations must be included. 
In general, each of these equations will have coefficients of all of the interfacial displacements, 

w k ( k = 0,1 N ), e 0 , and y 0 as well as pseudo force terms due to the nonlinearities and a 

"force" due to a change in temperature. The integral force equation will be of the form 

N 

k w k + ^N+I 6 © + ®n+2?o = F X + f nl + Ft [2.60] 

k=0 

and the integral torque equation will be 

N 

X/W + 'Wo + 'Wo = T x + fJu + Ft- [2-61] 

k=Q 

For example, consider an arbitrary 3 layer tube. Assembly of the global stiffness matrix 
yields: 
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* 





* 

p -1 

K 11 

K 12 

0 

0 

K 13 

kJ 4 

w 0 

*21 

(K 22 + ^ 1 l) 

i /2 

l '12 

0 

(kJ 3 + k? 3 ) (k 2 4 + k?„) 

W-| 

0 

y/i- 

*■21 

(Km + K?,) k ? 2 

(K 2 3 + K 13 ) (K 2 4 + K u ) 

w 2 

0 

0 

K 3 

*21 

K 22 

K 23 

K 24 

w 3 

% 

«>i 

0 2 

*>3 

^N +1 

^N +2 

So 

*0 

% 


^3 

*N + 1 


Vo 

m 


Pi 


Ft" 


pi— 

' NL 

0 


Ft + + Fj“ 


p1+ . p2— 

' NL "h P NL 

0 

-Po 


Ft + + Fj“ 

F? + 

+ 

p2+ , p3— 
P NL+ P NL 

p3+ 

P NL 

F x 


Ft 


f f 

' NL 

Tx 

. 


Ft 


Fnl 


or in short-hand notation 


[M]{W} = {F + F nl +F t } 


where 


[M] = global stiffness matrix 


{W} = {w 0 w, w 2 . . . w N e Q y 0 } T 
{F} = {P,0 ...0 -P 0 F x T x } t 
(F nl ) = global nonlinear force vector 
{F t } = global force vector due to temperature changes 


[2.62] 


[2.63] 
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This system of equations consists of N + 3 equations that must be solved for the radial dis- 
placements at the interfaces, w k , the axial strain, e 0 , and the rotation per unit length, y 0 . For 
a large number of layers, the local/global stiffness matrix scheme offers a nearly 50% re- 
duction in the number of equations to be solved. 


2.4 Endochronic Theory 


The Endochronic Theory is a nonlinear constitutive model based on irreversible 
thermodynamics. It uses the concept of internal variables and is centered around a defor- 
mation scale ( or intrinsic time scale ) which is a material property. The deformation scale 
can either be dependent or independent of time. For this investigation, it will be assumed that 
the deformation scale is independent of time. This theory was developed by Valanis 17 18 to 
describe such phenomena as cross-hardening in tension due to torsion and formation of 
hysteresis loops during loading and unloading cycles in isotropic materials. This theory was 
later extended to include anisotropic materials by Pindera and Herakovich 18 . A set of 
constitutive relations to model the in-plane behavior of transversely isotropic materials was 
obtained as a special case from the three dimensional formulation in order to correlate the 
theory with experimental data. The three dimensional relations developed by Pindera and 
Herakovich were subsequently incorporated into a three dimensional finite element model by 
Mathison, et a!., 20 to examine the behavior of compression-loaded laminates with a hole. 

The equations used by Mathison, et al M are ( without change ): 

e l = Atjo-j ± ! I ay | ni [2.64a] 
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[2.64b] 


*2 = + 


*2 

B? 2 (z-z') n * 

0 


d <*2 

dz' 


dz' 


s 3 — + 


r 

vn 


. K 3(To 

B? 2 (z-z ') n! -^-d2' 


[2.64c] 


*1 


5cr> 


«4 = A 4J®)+ 2 | B£ s (Z — Z')" 2 dZ' 


[2.64d] 


I 


e 5 = Asj^j + | ®66( z “ z 0 n8 q z > dz * 


H 


= A 6 jff]+ I 
Jq 


Bee(z — z') n * ~^~ir dz' 


[2.64e] 


[2.64f] 


where B 11t B&, Bfo, n 1t n 2 , and n e are nonlinear material constants, A y are initial compliances 
of the material, t * is the variable of integration and 


dz = • v /S|jd<X|derj 


and Sy is a fourth order positive definite tensor whose components are material properties. 
Bf 2 and B& are constants pertaining to the dissipative type nonlinear response. B„ , on the 
other hand, describes the reversible nonlinear behavior in the fiber direction. 
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Recalling from the previous section, it was found that T xr = r flr = 0 throughout the thick- 
ness of the tube. Transforming these stresses to the principal material coordinates system, 
it can be shown that <r 4 = 0 and cr 5 = 0 throughout the thickness of the tube also. Extracting the 
nonlinear portion of the strain from equations [2.64d] and [2.64eJ 


NL 

e 4 


if 


a " 4 


dz ' 


dz' 


NL 

e 5 


= j ^6 6 (Z — Z')" 8 


d<r 5 
dz ' 


dz' 


Since <r 4 — <r s — 0, 




therefore 


which is the result that was used in the previous section. Also, because <r 4 = <r 6 = 0, the elastic 
portion of the strains in Equations [2.64d] and [2.64e] are zero, therefore, the total strains, 
e 4 and £ 5 are zero. Restating the endochronic constitutive equations 


e 1 = A 1 j<Tj±B 11 |<r 1 i ni 


[2.65a] 


e 2 = ^2j a l + 


f 


B? 2 (z-z') ni “-dz' 


[2.65b] 
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c 3 = A 3j crj + 1 Bj 2 (z-z')" ! 4?rdz' 

[2.65c] 

J 0 


«4 = 0 

[2.65d] 

o 

(1 

[2.65e] 


E&(z-z') n '-|p-dz' [2.650 

Because the deformation scale, z, is a function of the loading history, which is, in general, 
complex, equations [2.65a-f] cannot be integrated exactly. Therefore, for these equations to 
be used in an numerical model, they must be integrated incrementally. The incremental form 
of the expression for e B was derived by Mathison.et al., 20 as follows: 




Zn 



k=1 
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— V Sjj AcTj A<Tj k 


( k not summed ) 


N = number of increments 

The incremental forms for e 2 and e 3 are similar to equation {2.661. 


2.5 Failure 


For a complete investigation of the effect of combined thermomechanical loading on the 
global response of composite tubes, the prediction of failure within the tube using an existing 
failure criteria is desired. Three failure theories were evaluated for incorporation into the 
analytical model: Maximum Stress Theory, Tsai-Hill Theory, and Tsai-Wu Theory. Because 
consideration for interaction between the different stress components and for the difference 
in magnitudes of ultimate stresses in tension and compression are desired, the Tsai-Wu The- 
ory was selected for use in the analytical model. 

According to Tsai and Wu 21 , failure will occur when 

FjO’j + Fjj<7j(7j = 1 i,j = 1,2,3,. ...6 

where contracted notation has been used. F, and F fj are second and fourth order tensors re- 
spectively that are functions of the ultimate strengths of the material: 



F4 = F 5 = F6 = 0 
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F44 = F55 = 0 


where X T and X c are the tensile and compressive strengths in the fiber direction, Y T and Y c 
are the tensile and compressive strengths in the transverse direction, and S is the in-plane 
shear strength. For the purpose of this investigation, the coupling terms, F 12) F 13 , and F 23 , will 
be neglected. This failure theory employs the stresses at a given point in the tube to predict 
failure, therefore, the prediction of failure using this theory in the analysis of tubes will only 
yield individual ply failures rather that the total failure of the tube. 


2.6 Numerical Procedure 


2.6.1 Approximation of « t ML (r) 


As stated in an earlier section, the nonlinear strains, ef L (r), are approximated by a fifth 
order polynomial, 

£j NL (r) = a, + b,r + c,r 2 + d,r 3 + e/ + f t r 5 
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in each layer. The unknown coefficients will be evaluated using the Least Squares Method 
(see Appendix B). This method requires that strains at six points within each layer must be 
known. Therefore each layer must be divided up into five or more sublayers of the same fiber 
orientation. This facilitates the calculation of £| NL (r) at each of the interfaces of the sublayers 
and, hence, permits the calculations of the unknown coefficients of the polynomial. 


2.6.2 Iterative Technique 


Many techniques for solving a system of nonlinear equations are currently in use includ- 
ing incremental procedures, iterative procedures, and step-iterative or mixed procedures 20 . 
For the solution of this problem, a step-iterative scheme was chosen. This scheme is simply 
a combination of an incremental procedure and a functional iterative procedure. 

The functional iterative technique starts with an initial estimate for the displacements, 
{W}, and performs the following iteration procedure: 

[M]{W (r) ] - {F} + {F t } + {F NL ({W (r ~ 1) })} 

where r represents the iteration number. The iteration is continued until the solution of two 
successive iterations satisfies the convergence criteria. 

The incremental procedure breaks the total load applied into finite increments. For each 
load increment, the functional iteration technique is performed to arrive at a solution for each 
load increment. 

The global stiffness matrix, [M], contains only linear terms. Therefore it need only be 
decomposed once for the entire loading cycle. AH of the nonlinear terms are contained in the 
pseudo force vector {F Nl }, therefore only the right-hand side of the equation must be changed 
from one iteration to the next. 

A technique called scaled partial pivoting is used to solve the system of equations at each 
iteration. This technique is suitable for the case where the same stiffness matrix is used with 
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many force vectors as is the case in this situation. This procedure utilizes Gaussian elimi- 
nation in which rows are interchanged and scaled in order to pivot from a position that will 
cause the least numerical roundoff errors. The elimination procedure is performed once on 
the global stiffness matrix, [M]. The row Interchanges are stored, and at each iteration, they 
can be performed on the force vector. A solution for the displacements is then obtained from 
a back substitution procedure. This technique offers a significant savings in computing time 
over inversion of the stiffness matrix and subsequent multiplication with the force vector at 
each iteration 22 . 

The procedure that is followed in the step-iterative scheme is the constant stress formu- 
lation. For this formulation, the iterations are performed at essentially a constant stress level 
and the corresponding strains and displacements are determined for each load step. The it- 
erative procedure followed for the constant stress formulation is as follows: 

1. The system of equations is solved for w, e 0 , and y 0 . 

2. Convergence of solution is checked. 

3. The global strains, e x , e„ and y * are calculated from the radial displacements, 

e 0 and y 0 using equations [2.14a], [2.14b], [2.14c], and [2.14f]. 

4. The global strains are transformed into the principal material coordinates, and the stress 
in the principal material coordinates are determined from the constitutive equations, 
Equation [2.4a], using the nonlinear strains, eP L (r), from the previous iteration. 

5. The stress increments are determined from the difference of the current value of stress 
and the converged value from the previous iteration. 

6. The nonlinear strains are calculated from the stress increments. 

7. The nonlinear strains in each layer are approximated by fifth order polynomial. 
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8. The pseudo force vector, {F NL }, is calculated from the nonlinear strains 

This sequence is repeated until convergence of the solution is attained. A flow chart for 
the computer program is shown in Figure 2 


2.7 Program Verification 


In addition to performing the nonlinear tube analysis, the program developed is also ca- 
pable of linear analysis. To verify the results obtained using the local/global stiffness formu- 
lation, the results obtained from the linear portion of the program developed for this study 
were compared with results obtained from the computer program developed by Rousseau, 
et al 13 . 

An arbitrary unsymmetric, unbalanced laminate with a stacking sequence of 
[10/45/ -30/60/0/ —15/5/ —75/90/50/ -60/30/0/ — 25/82] T was selected. The material used for 
this analysis was AS4/3502 graphite/epoxy (using the material properties given by Mathison, 
et al., 20 ). A general loading including all types of loading was chosen. The applied loads are: 

F x = 10000 lb. 

T x = 1500 in-lb. 

Pi = 225 psi 
P 0 = 100 psi 
AT = — 100°F 

A comparison of the distribution of stresses through the thickness is shown in Figure 3 
through Figure 6. There are no differences in the results from the two models. 

To verily the nonlinear portion of the program, a pure shear state can be used as a test 

case 
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Oq- 5*0 

<T| = 0, i = 1,2, 3,4,5. 


For the pure shear state, Equation [2.64f] can be integrated exactly 


e 6 — Aej<?j + 



B 6 6(z-Z') ne 


dff 6 

dz' 


dz' 


where 


therefore 


dz = % /S i jd<T|daj 
— J d<Tg 


[2.67] 


dz ~ dz 


1 



Substituting for 


da t 

dz 


in Equation [2.64f] 


( % z 


e 6 ~ ^6j° r j + 


*/q 


■J$66 


86 6JffJ+ ^(n 6 + 1) 


[ 2 . 68 ] 


Integrating Equation [2.67], we have z = -Js^a 6 . Substituting for z into Equation [2.68] 


Analytical Model 


46 



[ 2 . 69 ] 
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The shear strains can now be determined for any shear stress. 

A state of pure shear can be attained In a tube by applying pure torsion to a tube in which 
all of the layers are aligned with the fibers in the axial direction ( 0° ). Comparing the analysis 
of a 0° tube in pure torsion with results from Equation [2.69], the analytical model can be 
verified for the nonlinear analysis of tubes. These results are shown in Figure 7. 
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3.0 Experimental Procedure 


3.1 Introduction 


The experimental program consists of a series of flat coupon tests followed by biaxial 
tests on composite tubes. The flat coupon tests are comprised of a set of off-axis tension and 
compression tests and losipescu shear tests on coupons cut from a flat, 12 ply unidirectional 
panel made from P75/934 graphite/epoxy. The results of these tests will be used to generate 
material parameters that will subsequently be used in the analytical model to predict the re- 
sponse of the composite tubes. The laminated tubes will then be tested in combined axial and 
torsional loading, and the results will be compared with predictions of the model. The extent 
of the test program on the laminated tubes was limited because fewer tubes were available 
than initially planned for. A test fixture was also developed for the application of combined 
loading. The description of the design of the test fixture will be included in this chapter. 
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3.2 Test Plan 


3.2.1 Material Characterization 

The test matrix for the tests performed on the flat coupons is shown in Table 1. The 
material used in this study will be assumed to be transversely isotropic so there are five 
elastic properties to be determined: E^, E 22 , v 12 , v 23 , and G 12 . The off-axis tension tests were 
performed to determine the in-plane material properties in tension and shear. Fiber orien- 
tations of 0° and 90° were used to determine the axial and transverse Young's moduli as well 
as the Poisson's ratio, v 12 . Fiber orientations of 10° and 45° were used to determine the axial 
shear modulus, G 12 . The 45° orientation provides the most accurate prediction of the shear 
modulus while the 10° orientation gives a measure of the effect of the Interaction of the 
i and <r 2 stress components on the shear response of the material. These tests were also 
used to determine the ultimate strengths of the material in tension. 

The off-axis compression tests were performed to characterize the in-plane material re- 
sponse in compression. Three fiber orientations were tested: 0°, 45°, and 90°. The ultimate 
strengths in compression were also determined. 

The losipescu tests were performed to characterize the in-plane shear response of the 
material under nearly pure shear conditions. The results from these tests will be compared 
with the 10° and 45° off-axis tests. Fiber orientations of 0° and 90° were used for the losipescu 
tests. 

To supplement the tests performed on the flat coupons, tests were performed on 
unidirectional tubes. The test matrix for these tubes is shown in Table 2. Tubes with fiber 
orientations of 0°, 45°, and 90° were tested as a comparison for the flat panel tests. 
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Table 1. Flat Coupon Test Matrix 


Test 

Fiber 

Number of 

Type 

Orientation 

Specimens 

Tension 

0° 

3 


10° 

3 


45° 

3 


90° 

3 

Compression 

0° 

3 


45° 

3 


90° 

3 

iosipescu 

0° 

3 


90° 

6 
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Table 2. Unidirectional Tube Test Matrix. 


Tube No. 

Type of Loading 

0-01 

Tension 

0-02 

* 

0-03 

Torsion 

45-01 

Torsion 

45-02 

* 

45-03 

Tension 

90-01 

** 

90-02 

Torsion, Tension *** 

90-03 

Torsion, Compression *** 


Received in damaged condition 
Failed in test fixture prior to testing 

Two tests accomplished by performing first test well below failure 
load then performing second test. 
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3.2.2 Laminated Tube Tests 


The laminated tube tests were performed on seven tubes with stacking sequences of 
[15/0/ ± 10/0/ — 15],. Four of the seven tubes had been thermally cycled between -150°F and 
+ 150°F for 3000 cycles prior to testing while three had not. The response to pure axial loading 
( tension or compression ) was found to be linear to failure in preliminary tests performed on 
laminated tubes. The pure torsion response, on the other hand, was nonlinear. Because the 
nonlinear effects are of primary interest, the emphasis on the tests on the laminated tubes 
will be in pure torsion and torsion dominated biaxial loading. The seven tubes were tested 
using four different loading types: 


Type I Loading Monotonic negative torsion to failure 

Type li Loading Cyclic positive to negative torsion. The magnitude of torque is increased 
with each successive cycle. 

Type ill Loading Combined Loading - tension/positive torsion. The tension and torque are 
increased proportionally to a specified magnitude then unloaded propor- 
tionally. Pure torsion is applied to the same magnitude then pure tension 
is applied. Finally the tube in unloaded proportionally. The ratio of torque 
to tension is 1.73:1. This procedure is repeated with increasing magni- 
tudes of loads. 


Type IV Loading Combined Loading - compression/negative torsion. The same loading 
sequence as Type III is used except it is loaded in compression and neg- 
ative torsion rather that tension and positive torsion. This procedure is 
also repeated with increasing magnitudes of loads. 
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Figure 8 shows the loading sequences used for each type. The test matrix for the laminated 
tubes is shown in Table 3. Table 4 through Table 6 show the magnitudes of the loads applied 
for Types II, III and IV. The designation for the tubes will be in the form LMXYY where 

• LM defines the stacking sequence as [15/0/ ± 10/0/ — 15] s 

• X = C represents a thermally cycled tube 

• X = B represents a baseline tube (no thermal cycling) 

• YY is a sequential numbering of the tubes 


3.3 Material 


The material system chosen for this study consists of a P75 graphite fiber made by Amoco 
preimpregnated in a matrix of Fiberite 934 epoxy. The P75 fiber is a high modulus fiber, and 
the 934 epoxy is a standard epoxy used in many aerospace applications. Prior to this study, 
only some preliminary material properties for the linear region were available so the material 
system was characterized for the in-plane materia! properties in both the linear and nonlinear 
regions before conducting the biaxial tests on the laminated tubes. 


3.3.1 Panel Configuration 


For the purpose of characterizing the material, flat panels were fabricated using a 12 ply 
unidirectional layup of the P75/934 graphite/epoxy. The 12 ply configuration was chosen to 
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Table 3. Laminated Tube Test Matrix 


Tube 

Condition 

Load Type 

LMB01 

Baseline 

IV 

LMB02 

Baseline 

II 

LMB04 

Baseline 

II 

LMC05 

Thermally Cycled 

II 

LMC06 

Thermally Cycled 

III 

LMC07 

Thermally Cycled 

IV 

LMC08 

Thermally Cycled 

1 
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- Type II Loading 


Tube 

LMB02 

LMB04 

LMC05 



5 



Table 5. Load Magnitudes - Type III Loading 


Tube 

LMC06 


Cycle # 


1 

2 

3 




2402 4160 

2546 4410 
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Table 6. Load Magnitudes > Type IV Loading 


Tube 

Cycle # 

Axial (lb) 

Torque (in-lb) 

LMB01 

1 

-640 

-1160 


2 

-815 

-1410 


3 

-965 

-1660 


4 

-1105 

-1910 


5 

-1250 

-2160 


6 

-1400 

-2410 


7 

-1535 

-2660 

LMC07 

1 

-640 

-1160 


2 

-815 

-1410 


3 

-965 

-1660 


4 

-1105 

-1910 


5 

-1250 

-2160 


6 

-1400 

-2410 


7 

-1535 

-2660 


8 

-1680 

-2910 
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match the thickness of the tubes. The average thickness for the panels was measured to be 
0.065". 

Micrographs of the panel's cross section were taken to measure the volume fraction of 
the panels as well as to examine the distribution of the fibers in the matrix. A micrograph is 
shown in Figure 9. The fibers are evenly distributed through the thickness, and the ply 
boundaries are not distinguishable except for a few locations. The volume fraction of the 
panels were determined from the micrographs using an optical technique 23 . This technique 
yielded 65% for the fiber volume fraction. 


3.3.2 Tube Configuration 

The tubes used in this study are those that have been proposed by NASA for the truss 
structure of the Space Station. These tubes have been constructed from P75/934 
graphite/epoxy with a stacking sequence of [15/0/ ± 10/0/ —15],. The nominal inside diameter 
of the tubes is two inches with a twelve ply construction giving a wall thickness of approxi- 
mately 0.060". For the purposes of this test program the tubes have been cut to a length of 
ten inches. In addition to the [15/0/ ± 10/0/ — 15] s configuration, which will hereafter be re- 
ferred to as the "laminated tubes", unidirectional tubes were also constructed in fiber orien- 
tations of 0°, 45°, and 90° as mentioned earlier. 

The tubes that were employed in this experimental program were obtained from two 
sources: Boeing Aerospace and Morton Thiokol. The tubes constructed by Boeing were pri- 
marily used for preliminary testing to evaluate the newly developed test fixture to be dis- 
cussed in Section 3.4.3.4, while the tubes made by Morton Thiokol were used to obtain the 
results presented in Section 4.3. The manufacturing techniques employed by the two suppli- 
ers were somewhat different. The tubes made by Boeing Aerospace were fabricated using a 
unidirectional tape layup on a male mandrel. In addition, a two mil layer of aluminum was 
co-cured onto the inner surface of the tube. Morton Thiokol, on the other hand, fabricated the 
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Figure 9. Micrograph of Panel 
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tubes using a female mandrel in which a tape layup was cured. These differences in manu- 
facturing techniques led to differences in the dimensional repeatability of the tubes. The 
Boeing tubes possessed consistent inside diameters from one tube to the next due to the use 
of an internal male mandrel. The Morton Thiokol tubes had very consistent outside diameter 
owing to the female mandrel used. Because the outside diameter on the Boeing tubes and 
the inside diameter on the Morton Thiokol tubes were not controlled by a tool surface, the 
variation in the dimensions of the respective diameters varied considerably in some cases. 

In addition to the laminated tubes, Morton Thiokol manufactured the unidirectional tubes. 
These tubes were made with a 12 ply layup using the same manufacturing technique and 
tooling as was used to produce laminated tubes with roughly the same dimensions. 

Micrographs of the cross sections of the 0° and laminated tubes were taken. These ap- 
pear in Figure 10 and Figure 11. The fiber volume fraction of the 0° tubes is significantly 
lower than that of the flat panels; the fibers are visibly much less dense than in the panel. In 
addition, the ply boundaries are much more pronounced. The fiber volume fraction for the 0° 
tubes was calculated to be 50% using the optical technique. The micrographs of the lami- 
nated tube show many of the same characteristics as the 0° tubes: pronounced ply bounda- 
ries and lower fiber volume fraction. It should also be noted that the plies toward the center 
of the laminate are compacted much less than the plies at the inner and outer surfaces. The 
overall fiber volume fraction was determined to be 51%. However, if the fiber volume fraction 
is determined within the layer only - not including the resin-rich regions at the interfaces - the 
volume fraction increases to 57%. The resin-rich regions at the ply interfaces, if considered 
as a separate ply, were determined to be approximately seven percent of the overall thickness 
of the laminate. With this knowledge, the "layers" of resin can be Incorporated into the ana- 
lytical model as separate layers. 

The fiber volume fractions of the unidirectional tubes and the flat panel were also deter- 
mined using the matrix digestion method according to the ASTM D-3171 test method. The 
volume fractions for the 0°, 45°, and 90° tubes were determined to be 53%, 52%, and 54.1% 
respectively 24 , and the fiber volume fraction of the panel was found to be 54.3%. The result 
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Figure 11. Micrograph of Laminated Tuba 
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for the 0° tube compares favorably with the result using the optical technique. The result for 
the panel, however, was significantly lower than was determined using the optical technique. 
According to the results using the matrix digestion technique, the fiber volume fractions of the 
0° tube and the panel are nearly the same, but a comparison the micrographs of the cross 
sections indicates that the fiber density in the panel is much higher than the 0° tube. There- 
fore the fiber volume fraction of 65% determined using the optical technique will be the value 
used in future analysis. 


3.4 Test Method 


3.4.1 Specimen Preparation 


Prior to testing, all specimens were ultrasonically inspected to check for large voids or 
flaws. The flat panels were inspected prior to cutting the test coupons so that voids or flaws 
in the panels could be avoided. No detectable flaws were found in any of the panels. The 
tubes were inspected to determine their relative condition prior to testing since flaws could 
not be avoided due to the limited number of tubes available. The laminated tubes showed 
signs of possible voids and delaminations. A typical C-scan of a laminated tube is shown in 
Figure 12. The light areas signify an area of a possible delamination or void while the dark 
areas denote an area of good bonding. The micrographs presented in the previous section 
do not, however, show signs of delaminations. Following the C-scanning procedure the 
specimens were placed in a vacuum chamber to remove the moisture absorbed during the 
C-scanning procedure. 

The off-axis tension coupons were cut on a diamond impregnated saw to nominal di- 
mensions of Vz" x 12" with fiber orientations of 0°, 10°, 45°, and 90°. Allowing for two inches 
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for grips at either end of the specimen, an aspect ratio of 16 was realized. At the center of 
each specimen a WK-00-120WR-350 Micro-Measurements stacked rosette was mounted with 
a WK-06-1 25AD-350 uniaxial strain gage mounted back-to-back to allow for bending correction. 
Three specimens of each fiber orientation were prepared. The off-axis compression speci- 
mens were cut to nominal dimensions of 1" x 1 Vz", and the edges were ground for smoothness 
and parallelism. A FRA-2-11 2mm Texas Measurements stacked rosette was mounted at the 
center of each specimen, and a FLA-2-11 2mm uniaxiai strain gage was mounted back-to-back. 
The losipescu specimens were cut and ground to the dimensions shown in Figure 13. A 
FRA-2-11 2mm stacked rosette was mounted at the center of the test section of each speci- 
men. 


3.4.2 Test Equipment 


Two testing machines were used for the experimental work in this investigation. For the 
off-axis tension and compression tests and the losipescu tests, a screw driven displacement 
controlled UTS Machine was utilized. The tube tests were performed on an Instron biaxial, 
servo-hydraulic testing machine. This machine was capable of performing tests in either 
displacement or load control; load control was used for the tube tests in this investigation. 

Computer controlled data acquisition systems were used with both the UTS and the 
Instron test machines. MATPAC2, a software package developed by Hidde, Beuth and 
Herakovich 25 was used in conjunction with the UTS machine. This system, in addition to ac- 
quiring the reading from the strain gages and load cell at three to four times per second, 
controlled the machine strain rate for the test. Following the test, MATPAC2 could correct for 
strain gage misalignment, initial curvature and bending, transverse sensitivity and actual fiber 
orientation. This software package was utilized on an IBM-XT Personal Computer interfaced 
with a Data Translation DT2805/5716 A/D - D/A board and a Vishay 2100 signal 
conditioner/amplifier. 


Experimental Procedure 


67 




Figure 13. losipescu Specimen Geometry 
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Data from the tests performed on the Instron test machine was acquired on a variation 
of the MATPAC2 software system entitled MATPACO. This system is similar to MATPAC2 
except that the Instron test machine was controlled independently of MATPACO. MATPACO 
was installed on an IBM-AT Personal Computer that was interfaced with an Orion Data Ac- 
quisition System signal conditioner. This system acquired data at a rate of two sweeps per 
second. 

3.4.3 Test Fixtures 

3.4.3.1 Off-Axis Tension 

The test fixture employed for the off-axis tension tests was a rotating end grip type fixture. 
This fixtures reduces the effect of the end constraint that is caused by preventing rotation in 
the grips of the fixture. For a detailed description of the test fixture and the test methods see 
Pindera and Herakovich 26 . 

3.4.3.2 Off-Axis Compression 

The fixture used for the compression tests employs an end loaded coupon supported by 
four circular pins to prevent out-of-plane displacements of the unloaded edges. A more de- 
tailed description of the fixture is given by Gurdal and Starbuck 27 and Pindera, et al. 28 
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3A.3.3 losipescu 


The test fixture used for the losipescu tests is a modified version of the original model 
developed by Walrath and Adams 29 . 

3. 4. 3. 4 Tube Test Fixture 

In order for composite tubes to be loaded in combined tension (or compression) and 
torsion, a suitable fixture for introducing the loads was designed. The problem of gripping 
composite tubes has previously been examined by Highton and Soden 30 and Toombes, 
Swanson, Cairnes 31 , and Guess and Haizlip 32 . Highton and Soden, working with a filament 
wound glass/epoxy tube, applied a tapered reinforcement consisting of additional layers of 
glass/epoxy and cast resin. A tapered collar was assembled around the reinforced end and 
bolted to flange at the end of the tube. Due to the taper, a radial stress was applied to the tube 
when bolted to the flange. This grip was designed for a specimen loaded in axial tension or 

j 

compression and subjected to internal pressure. 

Toombes, Swanson, and Cairnes 31 designed a grip similar to that of Highton and Soden 
for a tubular specimen to be loaded in axial compression and torsion. This design leaves no 
means for loading in axial tension. 

Guess and Haizlip 32 designed two different end grip configurations. The first consisted 
of a two part grip that is adhesively bonded to the inner and outer surfaces of the tube. Be- 
cause the adhesive did not provide sufficient maximum load a second design was considered. 
The second design consisted of an internal plug and an external split collar. A glass/epoxy 
fabric is overwrapped on the outside of the composite tube and cured onto the specimen. 
Threads are then machined onto the overwrap to mate with matching threads on the collar. 
The internal plug prevents collapse of the tube when the collar is clamped on the tube. This 
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design provides easy assembly and disassembly of the grips on the tube, but it requires that 
a glass/epoxy overwrap be applied to the tube and threads machined into the overwrap. 

With these designs in mind, some goals for a fixture design were established. First, the 
tubes should not require any major modifications ( such as was described by Highton and 
Soden and Guess and Haizlip ). Second, the tubes should be easy to change between tests 
so that many tubes can be tested in a short period of time. Third, the fixture should not 
damage the tube during removal from the fixture so that more than one test can be performed 
on the same tube. Finally, the fixture must be capable of introducing axial tension or com- 
pression in combination with torsion. From these goals, a set of basic design requirements 
was determined. 

The primary basis for the design of the fixture was that gripping the tube would be ac- 
complished using friction alone. In order to minimize the amount of pressure necessary to 
exert on the tube to grip it in the fixture, the load should be introduced on both the inner and 
outer surfaces of the tube and the amount of gripping surface should be sufficiently large. The 
fixture must allow for variations in the inner and outer diameters of the tube. 

The gripping fixture designed for this experimental program with the above criteria in 
mind mechanically grips the specimen by applying external pressure. The fixture consists of 
three basic components: the plug, the collar, and the base. All parts are made from carbon 
steel. The base consists of a circular plate and a protruding rod which can be gripped by the 
Instron hydraulic grips. The plug is bolted to the plate portion of the base using three 5/16" 
bolts and is then inserted into the end of the specimen. It is made in several different diam- 
eters so that a snug fit with the composite tube can be realized with the selection of the proper 
size plug. The collar, which consists of four identical parts, is then assembled around the 
specimen. It is also bolted to the base. Because both the collar and the plug are bolted to 
the base, the applied load is introduced to the specimen at both the inner and outer surfaces. 
This permits the gripping pressure to be reduced from what would be required if the load were 
introduced through only one surface. A schematic of the fixture is shown in Figure 14. 
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Figure 14. Tube Test Fixture Schematic 


The procedure of mounting the grip onto the specimen can be performed while the base 
at each end of the specimen is gripped in the Instron's hydraulic grips. This allows for auto- 
matic alignment of the tube with the axis of the testing machine and provides a stable platform 
in which the bolts in the fixture can be easily torqued without applying any load to the tube. 

A preliminary analysis was performed to determine the validity of the design outlined 
above. The length of contact between the tube and the collar and the plug was chosen to be 
1.75 inches. This corresponds to a contact area of 11.46 square inches on the inner and outer 
surfaces. Using the load capacity of the Instron testing machine (20,000 lb 
tension/compression and 10,000 in-lb torsion) as the maximum design load, the maximum 
shear stress required at the interface of the grip and the tube was determined to be 982 psi. 
A conservative value for the coefficient of static friction was assumed for the steel 
grip/composite tube interface to be ju = 0.20. From this, the magnitude of the radial stress 
generated by the grip required to deliver the necessary shear stress was determined to be 
4912 psi. 

Using the linear elastic portion of the analytical model developed in Chapter 2, the 
stresses induced into the tube by the grips alone were calculated. A worst case value of 9000 
psi for the gripping pressure was used. In addition to calculating the stresses due to the 
pressure exerted by the grips, calculations were made to determine the stresses when the 
plug inserted in the tube does not fit perfectly, i.e. there is a gap between the tube and the 
plug. It was found that for all but the case of a unidirectional 90° tube, the gripping pressure 
required to close a gap of 0.010" between the tube and the plug was small compared to that 
necessary to grip the specimen under the design loads. The 90° specimens, however, require 
grip pressures in excess of 40 ksi in order to close the gap due to the high stiffness in the 
circumferential direction. Therefore it is imperative the plugs match the inside diameter of the 
tubes perfectly ( or as near perfect as possible ) for the 90° tubes. Aside from the 90° 
unidirectional case, the stresses due to the gripping pressure were calculated to be well be- 
low the ultimate strengths of the material to be used ( P75/934 graphite/epoxy ). Because the 
main thrust of this work was not to develop an analytical tool for predicting the complex stress 
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states induced by this grip fixture, a more rigorous analytical method for this purpose was not 
pursued. The results of the performed analysis were determined to be sufficiently promising 
to further pursue this design. 

Once preliminary tests on the tubes manufactured by Boeing were initiated, it was noted 
that the tubes either slipped out of the grips or fractured in or near the grips. In order to al- 
leviate this, the surfaces of the plug and the collar that come into contact with the composite 
tubes were grit blasted to provide a rougher surface. This permitted gripping at lower levels 
of pressures exerted by the collar on the tube. The fracturing of the tubes in or near the grips 
was due to stress concentrations Induced by the grips. An attempt was made to alleviate this 
by using a torque wrench to tighten the grips onto the tube in a uniform pattern such that the 
collar is tightened evenly around the circumference of the tube. 

Several preliminary experiments were performed on the tubes fabricated by Boeing for 
the evaluation of the test fixture. Heavily instrumented tubes were tested in pure axial tension 
and pure torsion. For both axial and torsional loading, a six inch test section was found to be 
sufficient to minimize the end effects in order to obtain a uniform state of deformation midway 
between the grips. In pure tension, a maximum axial stress of 30 ksi was achieved with failure 
occurring in the fixture. This is 60% of the load for which first ply failure is predicted. Under 
pure torsion, failure of the tube occurred at an applied torque of over twice the predicted 
torque for first ply failure. These results indicate that this test fixture is suitable for use in 
determining the response and failure of tubes under torsional loading, but the testing of tubes 
under pure tension is limited to loads well below predicted failure. 
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4.0 Experimental Results 


4.1 Material Characterization 


The presentation of the experimental results will begin with a description of the initial 
elastic properties of the material followed by a discussion of the nonlinear properties. Be- 
cause the results of the characterization tests yielded consistent results, only representative 
test results will be presented in this chapter. A complete set of results is given in Appendix 
C. The data shown represents actual test data; no smoothing of the results has been per- 
formed. 


4.1.1 Tensile Properties 


Figure 15 shows the axial response of the four orientations of tensile coupons, and the 
Poisson's response for these tests is shown in Figure 16. The elastic properties and failure 
stresses obtained from these tests are summarized in Table 7. 
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Figure 15 . Axial Response of Tension Tests 
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Table 7. Summary of Average Tensile Properties 


Type 

Specimen 

E x (Msi) 

v *y 

G 12 *(Msi) 

cr^ LT (ksi) 

er<%) 

Panel 

1-00-01 

35.25 

0.340 


92.96 

0.225 


1-00-02 

35.25 

0.314 


103.69 

0.255 


1-00-03 

35.25 

0.340 


100.35 

0.256 


Average 0° 

35.25 

0.331 


98.76 

0.245 


1-10-01 

12.61 

0.450 

0.625 

42.31 

0.345 


1-10-02 

13.80 

0.357 

0.625 

43.10 

0.325 


1-10-03 

13.80 

0.366 

0.625 

32.73 

0.246 


Average 10° 

13.40 

0.391 

0.625 

39.38 

0.305 


3-45-01 

1.47 

0.309 

0.555 

7.75 

0.544 


3-45-02 

1.47 

0.331 

0.535 

5.91 

0.397 


3-45-04 

1.53 

0.295 

0.584 

7.30 

0.482 


Average 45° 

1.49 

0.312 

0.558 

6.99 

0.474 


3-90-01 

1.04 

0.016 


4.05 

0.388 


3-90-02 

1.04 

0.016 


3.76 

0.364 


3-90-03 

1.04 

0.016 


3.57 

0.346 


Average 90° 

1.04 

0.016 


3.80 

0.366 

Tube 

0-01 

36.38 


» - 

t 

t 


0-03 

46.67 



t 

t 


90-02 

1.10 

-•4 


t 

t 


90-03 

1.17 

-•4 


t 

t 


* Apparent 
| Tube not failed 

t Data too noisy for accurate measurement 
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Figure 17 and Figure 18 show a comparison of experimental resuits for Young's modulus 


and Poisson's ratio with predictions of transformation theory, equations [4.1] and [4.2], 


1 


Exx(^) ^11 


^ COS 4 <t> + ^ ^ 


2v 12 ^ _:_2 1 2 , , 1 _:_4 , 

— — I sin q> cos <p + — — sin <p 

G 12 E 11 J fc 22 


[4.1] 


v xy(^>) = Exx^) 


v 12 

Ell 


( sin 4 <£ + cos V) - ( -=r- + ^r— ) sin 2 <£ cos 2 <£ 

x X \ *=11 *-22 b 12 / 


[4.2] 


The stress state in a multilayered tube is, in general, three-dimensional, therefore, an 
out-of-plane Poisson's ratio, v 23 , is required for subsequent experimental/analytica! correlation 
to be presented in Chapter 5. Because only in-plane properties can be determined from the 
off-axis tests, a value for v 23 could not be determined experimentally. A value of 0.49 for v 23 
was assumed based on Datta et al 33 . 


4.1.2 Compressive Properties 


Figure 19 shows the results for the axial response of the compression tests performed 
on the three fiber orientations. The results for the Poisson's response from the compression 
tests are shown in Figure 20. A summary of the results for the material properties obtained 
from the compression tests is given in Table 8. Figure 21 and Figure 22 show a comparison 
of these results with transformation theory (Equations [4.1] and [4.2]). 

To insure that the specimens were failing in compression rather than buckling, the 0° arm 
of the rosette was plotted with the uniaxial gage that was mounted on the opposite face of the 
specimen. Each specimen exhibited similar behavior so only a representative curve will be 
presented. The front and back response of a 0° compression specimen is shown in 
Figure 23. A sudden drop in the load prior to where the two curves separate indicates that 
an initial compressive failure occurs before buckling initiates. Therefore the failure stresses 
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Table 8. Summary of Average Compressive Properties 


Type 

Specimen 

E x (Msi) 

% 

G 12 *{Msi) 

ar(ksi) 

*r(%) 

Panel 

C-00-01 

34.48 

0.314 


-37.27 

-0.115 


C-00-02 

36.30 

0.373 


-49.00 

-0.165 


C-00-03 

36.82 

0.363 


-47.20 

-0.160 


Average 0° 

35.87 

0.363 

-.- 

-44.49 

-0.147 


C-45-01 

1.49 

0.447 

0.483 

-12.50 

-0.664 


C-45-02 

1.70 

0.445 

0.600 

-18.84 

-1.350 


C-45-04 

1.70 

0.573 

0.486 

-24.84 

-1.760 


Average 45° 

1.63 

0.498 

0.523 

18.73 

1.258 


C-90-01 

1.20 

0.0064 


-28.09 

-2.584 


C-90-02 

1.20 

0.0072 


-27.50 

-2.411 


C-90-03 

1.20 

0.0100 


-27.50 

-2.295 


Average 90° 

1.20 

0.0079 


-27.70 

-2.43 

Tube 

0-03 

45.57 

«• — 

... 

t 

t 


90-03 

1.06 

-•-t 


t 

t 


* Apparent 
t Tube not failed 

t Data too noisy for accurate measurement 
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Figure 21. Comparison of Compressive Young's Modulus with Transformation Theory 
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observed in these tests are indeed representative of a compressive failure rather than a 
buckling failure. 


4.1.3 Shear Properties 


The shear modulus can be determined from 10° and 45° off-axis tension, 45° off-axis 
compression or the 0° and 90° losipescu tests. As described by Pindera et a!. 34 , the values 
of the initial shear moduli from the different tests must be corrected. The off-axis tension 
specimens are subjected to an end constraint effect because the ends of the specimens are 
rigidly clamped. The effect of the end constraint can be corrected by employing the Pagano- 
Halpin model as described by Pindera and Herakovich 28 . 

The losipescu specimens, although subjected to a nearly pure state of shear stress, have 
stress distributions across the test section that are not uniform. Because of this nonuniform- 
ity, the apparent values of G 12 calculated using the average applied shear stress will not be 
the true shear modulus. By carrying out a finite element analysis to determine the actual 
distribution of shear stress across the test section, a more accurate value for G 12 can be ob- 
tained by applying correction factors to the apparent data as pointed out by Pindera et al 34 . 
The results of these tests are shown in Table 9. Figure 24 shows results of the 10° and 45° 
off-axis and 0° and 90° losipescu tests before correction factors are applied. With the cor- 
rection factors incorporated, good agreement between these tests can be seen in Figure 25. 


4.1.4 Nonlinear Properties 


The constant B„ and the endochronic constants, Bf 2l B|e, n 2 , n 8 , and , S,j , are deter- 
mined by assuming that for uniaxial states of stress, the nonlinear portion of the strain, 6j NL , 
can be approximated by a power law 18 20 
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Figure 24. Shear Response Before Shear Correction Factors 
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Figure 25. Shear Response After Shear Correction Factors 
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Table 9. Initial Shear Properties 



t Pindera and Herakovich 26 

* Pindera et al 34 
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[4.3] 


g| NL = Cjfff*. (i = 1,2,. ..,6 i not summed ) 

This implies that the plot of In versus In o-, can be approximated by a straight line. Taking 
the natural logarithm of both sides of Equation [4.3] yields 

In £j NL = ln(Cj<7, ki ) 

= In C, + In af 1 ( i not summed ) [4.4] 

mi 

In £| — In Cj 4“ kj in <7j 

In the Cartesian coordinate system, Equation [4.4] can be represented as 

y = mx + b 


where 


w NL 
y = In £| 

x — In <T| 
m = k| == slope 
b = lnC ( = y-intercept 

Figure 26 shows the fiber direction response in tension and compression. The P75/934 
graphite/epoxy exhibits a stiffening response in tension while the response softens in com- 
pression. The stiffening behavior in tension is bilinear with an initial modulus of 35.25 Msi as 
shown in Table 7, and the final modulus is 43.70 Msi which is a 24% increase in stiffness. To 
more accurately model this material, both the stiffening and softening effects will be incorpo- 
rated. 

The response of the 90° uniaxial tests are shown in Figure 27. In tension, the 
<j 2 versus e 2 curve is linear to failure. In compression, the curve is linear for a significant 
portion of the response. For these reasons, the transverse response will be assumed to be 
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linear, therefore, the nonlinear endochronic constants associated with the transverse direction 
(B^, n 2 , S 22 ) will be set to zero. 

Figure 25 shows the corrected curves for the shear response from the off-axis tension 
and compression and the losipescu tests. While the 0° and 90° iosipescu test specimens are 
in a nearly pure state of shear, the 45° and 10° off-axis tensile specimens are not. Although 
there exist other significant components of stress along material principal directions in espe- 
cially the 10° off-axis test, the degree of nonlinearity of each of the curves is nearly the same. 
This observation leads to the conclusion that the nonlinearity in the shear component of stress 
is uncoupled from the other components of stress. In other words, the presence of 
cr 1t cr 2 , and a 3 does not contribute to the nonlinear response in shear at the material level. 

From the lamina tests discussed previously in this chapter, and, assuming the material 
is transversely isotropic, we conclude that the only nonlinear components of stress are 
<7, and <y 8 which are uncoupled. Because the nonlinear portion of the strain in the fiber di- 
rection is independent of the deformation scale, z, and the uncoupled nature of the nonline- 
arities, the only nonzero component of the Sy matrix is the S 68 component. With B£ 2 and n 2 
already determined to be zero, the only parameters required are B 11t Bg 8 , n 1t n 8 , and S 88 . 
However, since the nonlinear response in the fiber direction stiffens in tension and softens in 
compression, two sets of B n and n., will be necessary, i.e., for a, >0 

and for cr, < 0 

NL B c I ,n? 
e 1 -- | \ 

The nonlinear strain in the fiber direction is very well approximated by that of the power 
law, Equation [4.3] where C t = B}^ or B ft) and k t — nj( or nf) can be determined directly from 
a plot of In e^ L versus In a v nj(nf) is the slope and In Bj 1 ( In Bft) Is the y intercept. 
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The values of n 8t Bg 6 , and S 66 can be determined from the graph of In e£ L versus In <j 6 for 
the case of pure shear. For this stress state, the form of the nonlinear strains, e£ L is given in 
Equation [2.69] 


o — 

NL ®66(^66) 2 n B +1 

— "5i"i * 6 


[4.5] 


Using Equation [4.3], Equation [4.5] can be represented by the power law approximation 


NL 


e 6 




where 


. S!l 
r B 66( S 66) 2 
6 ~ n 6 + 1 

kg = n 6 + 1 

The value of n 6 can be determined directly from the slope of the In eg 11 - versus In <r e curve. 

Because both Bg 6 and S 66 must be determined from a pure shear test, they cannot be in- 
dependent quantities. Since Bg„ is simply a multiplicative constant, it will arbitrarily be set 
equal to 1. With Bg e known, S 66 can be determined from the y-intercept of the 
In eg 11 - versus In <r 6 curve. 


4.2 Summary of Material Properties 


With the determination of the endochronic parameters, the P75/934 graphite/epoxy has 
been completely characterized in both the linear and nonlinear regions. A summary of the 
material properties to be used in the analytical model is shown in Table 10. The coefficients 
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of thermal expansion presented In Table 10 were determined by Bowles 24 for P75/934 with a 
fiber volume fraction of 50%. 

The response and failure of the 0° specimens were significantly different in tension and 
compression. The 0° coupons exhibited a stiffening response in tension white the behavior in 
compression softened. In addition, the ultimate strength in compression, X c , was less than 
50% of the tensile ultimate strength, X T . The initial moduli of the tension and compression 
tests differed by only 2% so a common value of 35.25 Msi will be used for both tension and 
compression. 

The 90° specimens also exhibited a large difference in the failure properties in tension 
and compression; the ultimate strength in compression was over seven times the tensile ul- 
timate strength. The initial response for both tension and compression was found to be linear. 
Although the initial Young's modulus in tension and compression differed by 10%, a common 
value of 1.04 Msi will be used for simplicity in the analytical model. 

Under complex loading paths, unloading of the stresses can occur due to the sequence 
of loading. For an accurate analytical prediction of this kind of behavior, an appropriate 
scheme for modeling the unloading of the stresses is required. The two components of stress 
which exhibit nonlinear constitutive behavior for P75/934 are the shear component and the 
component in the fiber direction which were found to be uncoupled. Because the unloading 
response was not part of the material characterization performed in the experimental inves- 
tigation, assumptions will be made as to the nature of the two nonlinear components of stress. 

Because the nonlinear portion of the strain in the fiber direction was assumed to be in- 
dependent of the deformation scale, z, the response is reversible. The reversible behavior 
of graphite composites in the fiber direction was verified by Pindera and Herakovich 35 using 
cyclic tests. 

The nonlinearities in the shear component are dissipative in nature. For dissipative re- 
sponse, the stresses will be assumed to unload linearly while maintaining a constant nonlin- 
ear strain. 
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Table 10. Summary of Material Properties 


Elastic Properties 

E^ = 35.25 Msi 
E 2 2 - 1.04 Msi 
v 12 - 0.331 
v 2 3 = 0.49 
G 12 = 0.570 Msi 


Failure Properties 


X T 98.8 ksi 
X c = -44.5 ksi 
Y t = 3.8 ksi 

Y c = -2 7.7 ksi 
S = 5.85 ksi 


Thermal Properties 

= -0.584 x 10 _e /°F 
« 2 = 19.18x10 ~ 6 /°F 


Endochronic Parameters 


Bri 

= 0.7840292X 10" 12 

n} 

= 1.736261 

Bfi 

- 0.2541205X 10- 12 

nf 

= 2.056897 


= 1.00 

n 8 

= 2.859718 

S 6 e 

= 0.1221467X 10- 13 
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4.3 Laminated Tube Tests 


As was mentioned in the previous chapter, the laminated tubes were tested in pure 
torsion and in torsion dominated biaxial loading. For the purposes of this investigation, only 
the global shear response of the tubes will be presented. For both the analytical predictions 
shown in the next chapter and the experimental results, the shear response will be presented 
in graphs of t* vs. y* where is the average applied shear stress 

Vn 

T x0- j 

and J is the polar moment of intertia of the tube cross section. The shear strain, y*. is the 
strain on the surface of the tube either measured with a strain gage for the experimental re- 
sults or calculated using the analytical model. 

In the description of the experimental results that follows, the magnitude of the applied 
shear stress where the response begins to exhibit nonlinear behavior will be referred to as 
the "linear limit". This point is frequently called the proportional limit in the literature and 
textbooks. 


4.3.1 Type I Loading 


The response of tube LMC08 to pure monotonic negative torsion to failure is shown in 
Figure 28. The shear modulus of the tube in this test was G,* = 1.32 Msi. The stress at which 
the response deviated from linear response was -6.0 ksi. Failure occurred at -8.88 ksi al- 
though cracking was audible at approximately 8 ksi. 
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4.3.2 Type II Loading 


The shear responses for the cyclic torsion tests of tubes LMB02, LMB04, and LMC05 are 
shown in Figure 29 through Figure 31 respectively. Each figure contains all loading cycles for 
a given tube with each successive loading cycle offset from the previous one by 0.10% strain. 


4.3.2.1 Tube LMB02 


The results of the cyclic torsion test on tube LMB02 are shown in Figure 29. The shear 
modulus, Grf, is 1.42 Msi. The linear limit for this tube, as shown in the second and subse- 
quent cycles, is 4.0 ksi. 

On the final loading cycle, audible cracking in the tube was detected at an applied shear 
stress of approximately 9.5 ksi in positive torsion. On the subsequent load up in negative 
torsion, total failure of the tube occurred at a stress of 5.3 ksi which is approximately 60% of 
the stress attained in negative torsion in the previous load cycle. The unloading response in 
the final loading cycle exhibits noticeable hysteresis which may be a result of the cracking 
observed in the loading in positive torsion. 


43.2.2 Tube LMB04 

Figure 30 shows the shear response for tube LMB04. For both loading cycles, the shear 
modulus for the tube, G^, is 1.24 Msi. The linear limit, as observed in the first loading cycle, 
was 4.5 ksi in both positive and negative torsion. Failure occurred at the peak of the second 
cycle at an applied shear stress of 8.15 ksi. 

In the first loading cycle, the unloading from positive torsion occurred elastically. How- 
ever, in negative torsion, the response exhibited a hysteresis in the unloading. This type of 
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response could be due to either the onset of damage or a dissipative permanent strain in the 
material. On the subsequent loading in positive torsion, failure of the tube occurred at the 
peak of the loading cycle. 


4.3.2.3 Tube LMC05 

The shear response for Type II loading on tube LMC05 is shown in Figure 31. The Initial shear 
modulus for this tube was the same for positive and negative torsion in these tests: 
G* = 1.43 Msi. The response of this tube is linear up to an applied shear stress of 4.8 ksi at 
which point the response began to soften. Failure occurred in this tube in a similar manner 
to the failure in tube LMB02; damage initiated in the positive torsion load up followed by failure 
in negative torsion at an applied stress of -6.11 ksi which is less than had been achieved in 
positive torsion. This test did exceed the stress applied in negative torsion in the previous 
load cycle. 


4.3.3 Type II! Loading 


Tube LMC06 was tested using Type III loading; thirteen load cycles were successfully 
completed before failure occurred in the fourteenth cycle. The initial shear modulus of 1.27 
Msi, as well as the curvature of the response was the same for each increment of load, 
therefore only the two final load cycles from each test will be presented to show the response 
through the entire loading cycle and the characteristics of failure. The linear limit for both 
combined proportional loading and pure torsion was 3.5 ksi. 

Figure 32 shows the shear response for the thirteenth and fourteenth load cycles. The 
response due to combined proportional loading and unloading, as well as due to pure torsion 
and pure tension are indicated on the Figure. The initial shear response for both the combined 
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proportional loading and the pure torsion are the same. The response due to pure torsion 
softens more than the combined proportional loading. Tension applied in addition to the 
torsion stiffens the response. Following the application of tension, the applied loads on the 
tube are the same as had been applied in proportional loading. The shear strain from the two 
loading paths, however, differs by 0.01%. This difference in shear strain indicates that the 
response may be exhibiting path dependent behavior. 

The response of the final load cycle follows the same path in combined proportional 
loading until failure occurred at an applied shear stress of 9.5 ksi in combined proportional 
loading. 


4.3.4 Type IV Loading 


4.3.4. 1 Tube LMB01 


Tube LMB01 was tested in combined compression and negative torsion using the Type 
IV loading. Six load cycles were completed successfully before failure occurred in the seventh 
cycle. The shear response from the sixth and seventh load cycles is shown in Figure 33. The 
response for the seventh step is offset by 0.10% strain. The initial shear modulus for these 
tests was 1.37 Msi, and the linear limit for both pure torsion and combined proportional load- 
ing was -3.4 ksi. 

The response due to combined proportional loading and unloading are indicated in 
Figure 33 as are the responses due to pure torsion and pure compression. The addition of 
compression to the negative torsion tends to soften the response by a small amount. This is 
in contrast to Type III loading where the pure torsion response was softer than for the com- 
bined proportional loading. 
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Figure 32. Shear Response of Tube LMC06 - Type III Loading 
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For the final cycle the maximum loads were reached when loading proportionally in 
torque and compression. However, when loading in pure torsion, the tube failed at -5.6 ksi 
prior to reaching the torque that had been attained previously in combined proportional 
loading. The failure of the tube, in this case, is dependent on the loading path. 

4.2A.2 Tube LMC07 


Tube LMC07 was tested using the same procedures and load cycles as tube LMB01. 
Despite the thermal cycling performed on this tube, seven successful load cycles were per- 
formed with failure not occurring until the eighth cycle. Figure 34 shows the results from the 
final two load cycles. The initial shear modulus for these tests was 1.30 Msi, and the linear 
limit for both pure torsion and combined proportional loading was 4.2 ksi. 

The tendency for the pure torsion shear response to be stiffer than the combined loading 
as was detected in Tube LMB01 is not evident in this tube. The failure characteristic in this 
tube is the same as in LMB01: failure occurs when loading in pure torsion at -6.47 ksi fol- 
lowing the combined proportional loading and unloading where the torsion load had reached 
a higher level. 


4.3.5 Summary 


Tests of four different load types were performed on the laminated tubes. A summary 
of the results for these seven tests is shown in Table 11. 

For each of the seven tests performed on the laminated tubes, the initial shear modulus 
fell in the range between 1.42 Msi and 1.27 Msi. The experimentally observed variance in the 
modulus does not appear to be the result of any factor such as the effect of axial loading or 
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Table 11. 


Summary of Laminated Tube Tests 


Load 

Type 

Tube 

Gy# 

(Msi) 

Linear 

Limit 

(ksi) 

(ksi) 

yJTt ' 
<%) 

Loading 

at 

Failure} 

1 

LMC08 

1,32 

-6.0 

-8.88 

-0.7Q 

Pure 

Torsion(-) 

II 

LMB02 

1.38 

±4.0 

-5.35 

-0.36 

Pure 

Torsion(-) 

II 

LMB04 

1.36 

±4.5 

8.15 

0.88 

Pure 

Torsion(4) 

II 

LMC05 

1.42 

±4.8 

-6.11 

-0.42 

Pure 

Torsion(-) 

III 

LMC06 

1.27 

3.5 

9.5 

0.80 

Prqportiona! 

(4,4) 

IV 

LMB01 

1.37 

-3.4 

-5.6 

-0.41 

Pure 

Torsion(-) 

IV 

LMC07 

1.30 

-4.2 

-6.47 

-0.56 

Pure 

Torsion(-) 




t Corresponds to failure load rather than maximum load 

$ A single (4) or {-) represent? the sign of the pure torsion applied. 
The double { + ,40 represents the sign of the torque and axial force 
in the case of proportional loading. 
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differences between positive and negative torsion. This variance will therefore be attributed 
to experimental scatter. 

Two loading types, il and IV, were performed in which the test was duplicated for both 
baseline and thermally cycled tubes. For the Type II loading, two baseline tubes, LMB02 and 
LMB04, and one thermally cycled tube, LMC05, were tested. Comparing the test of the 
thermally cycled tube, LMC05, with the test on LMB04, it is evident that the thermally cycled 
tube progressed further in the loading cycle than the baseline tube. The second baseline tube, 
LMB02, successfully completed the load cycle in which both LMB02 and LMC05 failed, and 
failed on the subsequent load cycle. For the Type IV loading, the thermally cycled tube, 
LMC07, reached a higher applied shear stress at failure than the comparable baseline tube, 
LMB01. The results of these tests indicate that the effect of the thermal cycling on the tubes 
is not significant. 

The failure stress and the linear limit were dependent upon the type of loading. Under 
Type I loading, monotonic negative torsion to failure, the tube had a linear limit of -6.0 ksi and 
failed at an applied shear stress of -8.88 ksi. Four tubes of Types II and IV loading also failed 
in negative torsion, but they each had a loading history prior to failure. Under Type II loading, 
tubes LMB02 and LMC05 had linear limits of ±5.0 and ±4.8 ksi, respectively, and failure 
stresses of -5.35 and -6.11 ksi. In both cases of Type II loading, stresses of higher magnitude 
had been applied in positive torsion prior to failure. The effect of the loading history on these 
tubes was to reduce both the linear limits and the ultimate strength in negative torsion. 

Tubes tested in Type IV loading, LMB01 and LMC07, which was combined negative 
torsion/compression, also exhibited a reduction in the linear limits and the failure stress. Both 
tubes tested under this loading type failed under pure negative torsion after loading in com- 
bined proportional loading where the applied torsion had reached a higher magnitude. Prior 
to the load cycle in which the failure occurred, tube LMB01 had experienced six loading cycles 
while tube LMC07 experienced seven complete load cycles. The reduction in the failure stress 
appears to be a result of this loading history. 
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The linear limits experienced by tubes LMB01 and LMC07 under Type IV loading were 
over 30% less than the linear limit observed in tube LMC08 under Type I loading. Because 
the linear limits under Type IV loading were the same for each load cycle, the reduction from 
the linear limit under Type I loading appears to be the result of the addition of the pornpressive 
load rather than the loading history. 

Type III and IV loading were carried out using the same ratios of torsion to axial toad, and 
the magnitudes of the load increments were the same for each type. The only difference be- 
tween the two load types was the sign of the loading. Failure in the Type III loading occurred 
at 9.5 ksi which is approximately a 50% increase in magnitude of the failure in Type IV loading. 
In addition to the magnitude of the stress at failure, the loading at failure differed between the 
two types; Type III failed under combined proportional positive torsion/tension loading while 
the Type IV tubes failed under pure negative torsion. 

For the load levels tested, only the Type III loading exhibited a slight path dependent re- 
sponse. The path dependent phenomena did not become apparent in Type III until the mag- 
nitudes of the toads applied had surpassed those that had been attained under Type IV 
loading. 

When failure occurred In the tubes under any of the four loading types, a crack along the 
length of the tube would propagate along the length. At the free edge, the tube would exper 
rience an offset of the free edge in the axial direction. This offset is a result of residual 
stresses present in the tube due to the fabrication process. 
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5.0 Analytical / Experimental Correlation 


5.1 Introduction 

In this chapter, results from the laminated tube tests presented in Section 4.3 will be 
compared with predictions of the analytical model developed In Chapter 2. The material 
properties determined for P75/934 graphite/epoxy in Chapter 4 will be corrected for differ- 
ences in fiber volume fraction between the flat panels and the laminated tubes. These cor- 
rected properties will then be used as input for the analytical model predictions. The effect 
of residual stress on the initial response will be examined as well as the importance of the 
consideration of nonlinear material behavior at the ply level on the predictions of the tube 
response and failure. 
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5.2 Parameters Affecting Initial Response 


5.2.1 Material Property Correction 


As discussed in the Chapter 3, there is a discrepancy between the fiber volume fraction 
of the laminated tubes and the flat panels. The difference in the fiber volume fraction will re- 
sult in a difference in material properties. In order to accurately predict the response of the 
tube with the analytical model, an accurate set of material properties is required for the tubes. 

The causes for the discrepancy between the fiber volume fractions of the panel and the 
laminated tube are the lower density of fibers within the layers and the existence of resin-rich 
regions separating these layers. For accurate predictions from the analytical model, consid- 
eration for both causes will be made. The properties of the layers in the tube will be deter- 
mined with the help of a micromechanical model using the known material properties 
determined from the flat panel, and the resin-rich regions between the layer will be charac- 
terized for use in the analytical model. 

5.2.1. 1 Micromechanical Mode! 

The micromechanical model chosen for the purpose of correcting for the lower fiber vol- 
ume fraction is the model developed by Aboudi 3 ®. The fibers are assumed to be square and 
arranged in a uniform, square array. A representative cell is chosen that consists of four 
subcells; one subcell contains the fiber while the other three contain the matrix, see 
Figure 35. By requiring that the continuity of tractions and displacements at the interfaces of 
the subcells be satisfied, the average properties of the lamina can be determined from known 
fiber and matrix properties at a given fiber volume fraction. 
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Rather than determining the properties of the composite from known constituent proper- 
ties, the properties of the fiber will be determined from known matrix and composite proper- 
ties. The lamina properties are already known as determined in the previous chapter. The 
properties of the 934 epoxy matrix have been determined by Fox et al 37 . For the fiber volume 
fraction of the panels (65%), the fiber properties can be determined. The results of this pro- 
cedure are: 

E a = 53.90 Msi 
E t = 1.24 Msi 
v A = 0.343 
v-p = 0.440 
G A = 1.120 Msi 

where the subscript A represents the axial direction and T the transverse direction. The re- 
sults given above have been independently verified by Professor Aboudi 38 . 

From the known fiber and matrix properties, the lamina properties can now be calculated 
for a fiber volume fraction of the layers in the tube ( 57%) from the Aboudi Model. The re- 
duced lamina properties of the P75/934 graphite/epoxy in the tubes are: 

E^ = 30.99 Msi (-12.1%) 

E 22 = 1.003 Msi (-16.4%) 
v 12 = 0.352 ( 4*6.3%) 

v 23 = 0.499 ( +1.8%) 

G 12 = 0.515 Msi ( -9.6%) 

where the percentage changes from the properties determined from the material character- 
ization tests are also given. The nonlinear behavior of the material will be assumed to be 
unaffected by the change in fiber volume fraction. 
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5. 2.1. 2 Modeling the Tube Cross Section 

In Section 3.4, the configuration of the laminated tubes was discussed, ft was determined 
that the resin-rich regions at the ply interfaces and inner and outer surfaces, if considered as 
separate plies, accounted for approximately 7% of the overall thickness of the tube. 

So that the thickness of each resin layer and graphite/epoxy layer does not have to be 
measured from the micrographs for each tube for Incorporation Into the analytical model, se- 
veral assumptions will be made. Since all seven tubes that were tested were cut from a single 
ten foot long tube, the fiber volume fraction of the graphite/epoxy layers and the percentage 
of resin layers will be assumed to be the same for each tube. For simplicity, we will assume 
that all of the graphite/epoxy layers are the same thickness. The same assumption will be 
made for the resin layers. 

Using these assumption, the procedure for determining the thickness of each of the plies 
is as follows. The overall thickness of the tube is measured in four locations around the cir- 
cumference of the tube and averaged. The total thickness of the graphite/epoxy layers Is 
determined by multiplying the average overall thickness by 0.93; this is In turn divided by the 
number of layers (12) to obtain the individual ply thickness. Similarly, the individual ply 
thickness for the resin layers is obtained by multiplying the average overall tube thickness by 
0.07 and dividing by the number of resin layers (13). 

Figure 36 shows a comparison between the pure cyclic torsion experiment (Type II) on 
tube LMB02 and the predictions of the analytical model using uncorrected material properties 
(as determined in Chapter 4) and corrected properties and cross section in the absence of 
residual stress. The initial shear modulus for the experimental curve is 1.38 Msi, and the 
shear modulus for the analytical model using corrected and uncorrected material properties 
are 1.50 Msi and 1.32 Msi. The correction of the material properties and consideration for the 
actual cross section does more accurately predict the initial shear response of the exper- 
iment. 
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Figure 36. Comparison of Corrected to Uncorrected Properties 
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5.2.2 Residual Stress 
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As pointed out in the previous chapter, there is evidence of residual thermal stresses 
being present in the laminated tubes that were tested. For the purpose of predicting the 
magnitudes of the residual stresses on the basis of the outlined analytical model, a stress-free 
temperature of 350°F will be used. This temperature is the cure temperature for P75/U34 
graphite/epoxy, and, although it may be a conservative (high) estimate for the stress-free 
temperature, a more accurate estimate is presently not available. All of the experiments 
performed on the laminated tubes were carried out at room temperature so the temperature 
difference used to calculate the residual Stress is -275°F. 

According to linear analysis, the inclusion of residual stresses does rtot affect the initial 
response due to a mechanically applied load. Nonlinear analysis, on the other hand, does 
predict a difference between the initial modulus With fesiduai Stress and without residual 
stress. Figure 37 shows a comparison between pure cyclic torsion on tube LM802 and non- 
linear analysis with and without residual stresses. The incorporation of the residual stresses 
decreases the initial modulus of the shear response predicted by the analytical model from 
1.32 Msi to 1.20 Msi. 

The reason for the difference in the initial response predicted by linear and nonlinear 
analyses is the state of stress in the tube due to the temperature change. Figure 38 and 
Figure 39 show the distribution of the and c 8 stress components through the thickness of 
the tube. The solid lines represent the residual state of stress for AT ~-275°F, and the dashed 
lines represent the state of stress resulting from the subsequent application of the first load 
increment in torsion. It is evident that the external +15° layers are the most highly stressed 
layers with axial stress of -12.5 ksi and shear stresses of -1.7 ksl following the application of 
torsion. The result of the presence of residual stresses in these layers Is that the nonlinear 
region of the material is reached as soon as the torsion Is applied. 
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5.3 Predictions of the Analytical Model 

Nonlinearity in the response of the laminated tubes is evident from the results of the ex- 
periments presented in Chapter 4. The prediction of the nonlinear behavior is an important 
objective for the development of the analytical model. 


5.3.1 Stress-Strain Response 


5. 3.1.1 Contributions to Global Nonlinearity 

The nonlinear stress-strain response of the tube, prior to the occurrence of localized 
failure in the individual layers, is a result of the uncoupled nonlinear response in the axial 
direction - stiffening in tension and softening in compression - and to shear stress. Because 
the nonlinear stress components are uncoupled, the contribution of each nonlinear component 
to the global nonlinear shear response can be investigated. Figure 41 and Figure 40 show 
the global response of the tube under pure cyclic torsion loading (Type II loading) plotted with 
the separated effects of the axial and shear nonlinearities. Figure 40 shows the response 
without residual stress and Figure 41 shows the response with residual stress. Without resi- 
dual stresses, the initial moduli of the three curves are identical as demonstrated in the pre- 
vious section. In this case, the response in which the only nonlinearity considered is cr 1 
follows nearly the same path as the case for which both and cr 8 nonlinearities are consid- 
ered. The response due to the consideration of nonlinearities in <r 8 only is nearly linear. The 
same trends are evident for the case of residual stresses. This indicates that the global 
nonlinear response to torsion is due primarily to the nonlinear behavior In the axial direction, 
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and, since the response is softening, the nonlinearities are dominated by the nonlinear re- 
sponse in compression. 

5. 3. 1.2 Combined Loading Response 

Figure 42 shows a comparison of the predicted global stress-strain response to Type III 
loading using linear and nonlinear analysis. The linear curve predicts the same shear re- 
sponse for combined proportional loading and pure torsion. In other words, the linear model 
does not predict coupling between axial loading and shear strain. 

The nonlinear analysis predicts the same initial shear modulus for the combined pro- 
portional loading and for pure torsional loading which is softer than that predicted using linear 
analysis as demonstrated in Section 5.2.2. The predicted response due to combined propor- 
tional loading exhibits a higher degree of softening than for pure torsion. In addition, coupling 
between axial tension and shear strain is present. 

Figure 43 shows the linear and nonlinear analytical predictions of the global stress-strain 
response for Type IV loading. As with the predictions for Type ill loading, the linear model is 
incapable of predicting differences in response due to the different loading paths. For this type 
of loading, the nonlinear analysis predicts the response to pure torsional loading to be more 
nonlinear than the response due to combined proportional loading. This is in contrast to the 
predictions for Type ill loading. Coupling between axial loading and shear strain is also pre- 
dicted for this type of loading. 


5.3.2 Failure 

The nonlinearity observed in the global response of the tubes tested can be the result 
of material nonlinearities as well as damage evolution in the form of individual ply failures. 
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Figure 40. Contributions to Global Nonlinearity - Without Residual Stress 
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Figure 41. Contributions to Global Nonlinearity - With Residual Stress 
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Figure 43. Comparison of Linear and Nonlinear Analysis - Type IV Loading 
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In order to properly evaluate the underlying causes of the observed nonlinear response, it is 
necessary to predict the magnitude of the applied loads at which first ply failure occurs. Using 
the Tsai-Wu failure criterion described in Chapter 2, a failure envelope for the laminated tubes 
has been determined. The failure envelope will describe the state of applied average shear 
stress, t *9 , and average axial stress, <r x , at which first ply failure in the laminated tube will 
occur. 

Figure 44 shows the failure envelopes for the linear and nonlinear cases where residual 
stresses have not been included. These surfaces were determined using three points in each 
quadrant. Each surface in this figure is symmetric about r^O, but, due to the large differ- 
ence between the tensile and compressive ultimate strengths in the fiber direction, the failure 
surface predicts failure for larger magnitudes of combined loading for which the axial com- 
ponent is positive than for which the axial component is negative. The linear and nonlinear 
failure surfaces predict similar loading state at failure for combined loading states in which 
a x is compressive, but significant differences exist between the linear and nonlinear failure 
surfaces for combined loading states with positive o*. in addition to predicting lower failure 
loads for positive a x , the curvature of the nonlinear failure surface is concave. Palmer, Maier, 
and Drucker 39 discussed the issue of concavity in reference to yield surfaces. They found that 
for a composite for which one of the constituents is elastic-perfectly plastic and the other is 
nonlinearly elastic of the stiffening type, the yield surface in concave. They concluded that the 
concavity was a result of the stiffening behavior of the one constituent. Mathematically, the 
descriptions of both the yield surfaces and failure surfaces involve quadratic functions of the 
stress components whose constant values represent convex surfaces in the appropriate 
stress space. Depending upon the definitions of yield and failure that are employed, the yield 
and failure surfaces can be coincident. The concavity of the failure surface shown in 
Figure 44, then, is due to the stiffening behavior in the fiber direction. 

Figure 45 shows the results of adding the effect of residual stresses for the linear and 
nonlinear cases. The residual stresses have the effect of reducing the failure loads for ail 
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proportions of loads, and, interestingly, the concave characteristics disappear in the presence 
of residual stresses. 


5.4 Correlation of Results 


In this section the predictions of the analytical model for Types II, 111, and IV loading will 
be compared with the experimental results. Corrected material properties as determined in 
Section 5.2.1 will be used, and the effects of residual stresses will be included. 


5.4.1 Type II Loading 


A comparison between analysis and experimental results for the test on tube LMB02 has 
been shown earlier in this chapter in Figure 37 where the stress-strain response of the third 
of four load cycles is shown. As noted previously, the prediction of the Initial modulus ( in- 
cluding residual stress ) is 1.20 Msi, and the shear modulus determined from the experiment 
is 1.38 Msi. In terms of nonlinearity, the experimental results exhibit a slightly larger amount 
of softening than is predicted by analysis. The linear limit from the experiment, as shown in 
Table 11 is 4.0 ksi; the analytical prediction for the linear limit is 4.6 ksi. 

The nonlinear analysis predicts first ply failure in the -15° layers at 4.17 ksi in posi- 
tive torsion. In negative torsion, first ply failure is predicted in the +15° layers at t* = -4.23 
ksi. The onset of the experimentally observed nonlinearity appears to correlate well with the 
analytical prediction for first ply failure. Tube LMB02 failed in negative torsion at r * = -5.35 
ksi subsequent to loading in positive torsion to higher magnitudes of applied shear stress. 
This failure stress is also well below that achieved in monotonic negative torsion, -8.88 ksi. 
Tube LMC05 exhibited similar characteristics of failure. The reason for failure at an appar- 
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Figure 45. Failure Envelope Based on Linear and Nonlinear Analyses - With Residual Stress 
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ently low applied stress is that In previous load cycles, stress at first ply failure had been ex- 
ceeded for both positive and negative torsion. Therefore both the +15° layers on the surface 
of the tube and the -15° layers on the interior of the laminate have been failed, reducing the 
load bearing capability under cyclic loading in, say, negative torsion only. 


5.4.2 Type III Loading 


Figure 46 shows the comparison between analysis (shown previously in Figure 42) and 
the experimental results (shown previously in Figure 32) from the test performed on tube 
LMC06 under combined tension/positive torsion. The curves shown represent the results from 
the thirteenth of fourteen cycles. The prediction of the initial shear modulus matches the ex- 
perimentally determined shear modulus of 1.27 Msi. The linear limit determined from exper- 
iment is 3.5 ksi and 5.5 ksi determined from the analysis. The trends predicted by the 
analytical model correlate well with the experiment; the pure torsion response exhibits a 
higher degree of nonlinearity than the combined loading and coupling between the axial ten- 
sion and shear strain does exist. At peak loading of the experiment, a difference In shear 
strain of 0.01% was observed between the two loading paths. This path dependence is not 
predicted by the analytical model. 

Under proportional loading, nonlinear analysis predicts first ply failure at t* = 4.46 ksi. 
The prediction for the first ply failure under pure torsion is —4.17 ksi. These are indicated 
in Figure 46. For both paths the failures are predicted to initiate in the -15° layers in the center 
of the laminate. Because the experimentally observed linear limit, as well as the predictions 
for first ply failure, fall below the prediction for the linear limit, the onset of nonlinearity may 
be the result of first ply failure. Experimentally, the tube failed under combined proportional 
loading although the first ply failure is predicted at a lower stress under pure torsion. For 
stress levels less than the predicted first ply failure, the agreement between analysis and 
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Figure 46. Analytical/Experimental Correlation • Type III Loading 


Analytical / Experimental Correlation 


135 




experiment is good. At stresses higher than the predicted first ply failure, the experiment 
becomes more nonlinear than is predicted by the analysis. 


5.4.3 Type IV Loading 


Figure 47 shows a comparison between analysis and the results from the experiment 
performed on tube LMB01 under combined compression/negative torsion. The experimental 
results are from the sixth of seven load cycles. The initial shear modulus for this test, as 
shown in Table 11, is 1.37 Msi while the analysis predicts 1.22 Msi. For this test, the initial 
shear modulus is more closely predicted using the linear analysis. This indicates that the 
prediction of the initial modulus using nonlinear analysis in the absence of residual stress 
would more closely match the experimentally observed shear modulus. The curvature of the 
prediction, however, matches that from the experiment. The trends exhibited in the exper- 
iment, stiffer response in pure negative torsion which is softened with the application of pure 
compression, is predicted by the analytical model. The experimentally observed linear limit 
is -3.4 ksi. Using nonlinear analysis, -4.0 ksi is predicted for the linear limit. 

First ply failure is predicted at t* = -3.77 ksi under proportional loading and t* = -4.28 
ksi under pure torsion. These are indicated in Figure 47. For both the pure torsion and 
combined loading, the failure is predicted in the +15° layers which are on the outside and 
inside surfaces of the tube. The dominant stress in the failure prediction is the <r A which is 
compressive under these loading conditions. For both tubes tested under Type IV loading, 
failure of the tube occurred when loading in pure torsion, although the lowest prediction for 
first ply failure is for proportional loading. 
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Figure 47. Analytical/Experimental Correlation • Type IV Loading 
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5.5 Discussion 


In two of the tubes tested under Type II loading, failure occurred in negative torsion at 
magnitudes of stress lower that had been reached previously. This apparent premature fail- 
ure can be attributed to the ply failures occurring in both the -15° layers and the +15° layers 
due to the reversal of the torsion. 

Under Type III loading, there are differences in the nonlinearity of the analytical prediction 
and the experimental results. In addition, a difference of 0.01% strain was observed exper- 
imentally between the two loading paths, but the analysis did not predict it. These differences 
can be explained by an accumulation of damage in the tube. The presence of damage may 
enhance the nonlinearity of the material and induce a path dependence on the material re- 
sponse. If damage is occurring in the tube, the nonlinear response is no longer a pure ma- 
terial response but a combination of material and structural response. 

The above hypothesis is supported by the correlation between the predictions of first ply 
failure and the experimentally observed initiation of nonlinearity in the global shear response. 
Table 12 summarizes the results presented in Section 5.4. Because the initiation of damage 
can both the reduce the linear limit and enhance the nonlinearity, the consideration of damage 
is necessary for accurate prediction of nonlinear response. 

As stated previously, the difference between Type III and Type IV loading was only in the 
sign of the loads applied; the proportions of torsion and axial loading were identical. The 
applied shear stress at the failure of the tube was significantly less for Type IV than for Type 
III. The first ply failure predictions for the different loading types were within 15% of one an- 
other, and the mode of failure for both types was compression in the fiber direction. The dif- 
ference between the two loading types is in the location of the first ply failure prediction. For 
Type III loading, first ply failure is predicted in the -15° layers in the center of the laminate, 
and the first ply failure is predicted in the external +15° layers for Type IV loading. This result 
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Table 12. Comparison Between Linear Limits and First Ply Failure 


Load 

Type 

Tube 

Linear 

Limit 

(Experiment) 

(ksi) 

Linear 

Limit 

(Analysis) 

(ksi) 

FPF 

(Nonlinear) 

(ksi) 

FPF 

(Linear) 

(ksi) 

II 

LMC02 

±4.0 

±4.6 

4,17 

-4.23 

4.16 

-4.05 

III 

LMC06 

3.5 

5.5 

4.46(prop) 

4.17(tor) 

4.48(prop) 

4.16(tor) 

IV 

LMB01 

-3.4 

-4.0 

-3.77(prop) 

-4.23(tor) 

-3.78(prop) 

-4.05(tor) 
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indicates that a ply failure on the external surface of the type may be more critical to the in- 
tegrity of the tube than a failure in an internal layer. 

In Section 5.3, it was pointed out that the dominant component of stress in the global 
nonlinear response of the tube was the a, component. The cr t component of stress also is the 
dominant component in the predicted failure of the tubes. Both the softening nonlinearity and 
failure in the fiber direction are a result of the compressive properties of the material. The 
cyclic pure torsion tests of the Type il experiments presented in Chapter 4 exhibited a non- 
linear reversible unloading response up until the load cycle in which failure occurred. The 
experiments performed using Type III and Type IV loading also exhibited a reversible behavior 
with the successive loading and unloading cycles. These reversible characteristics of the 
experiments support the assumption that the response in the fiber direction, whether stiffening 
in tension or softening in compression, is elastic. Prior to the failure observed in Type II 
loading, the stress-strain response exhibited dissipative behavior. In view of the fact that the 
first ply failure had been exceeded for both positive and negative torsion, it appears that this 
dissipation is the result of damage that had accumulated in the tube. 
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6.0 Conclusions and Recommendations 


The results of this study indicate that the consideration of material nonlinearities has a 
significant effect on the predicted response and failure of tubes made from P75/934 
graphite/epoxy with a stacking sequence of [15/0/ ± 10/0/ —15],. The softening response of 
the individual plies in the fiber direction under compressive stress was found to be the domi- 
nant factor in the observed global nonlinear response of the tubes tested under the loading 
conditions employed in this investigation. The softening in shear was found to make very little 
contribution to the global nonlinear response. The stiffening response in the fiber direction 
due to tension results in the curvature of the failure surface becoming concave in the absence 
of residual stress. Nonlinear analysis was found to be capable of qualitatively predicting the 
experimentally observed differences in the shear response under pure torsion and combined 
proportional loading. The observed coupling between axial loading and shear strain was also 
predicted using the nonlinear analysis. However, the correlation between the initiation of 
global nonlinear response and the first ply failure predicted by the analytical model indicates 
that the observed nonlinear response is due to both material nonlinearity at the ply level and 
damage evolution due to successive ply cracking. 

The inclusion of residual stresses in the determination of the response and failure was 
found to be important. In the case of nonlinear analysis, residual stresses reduce the initial 
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shear modulus of the tube. The magnitudes of failure loads were also reduced by the residual 
stresses, and the concavity of the failure surface mentioned above was eliminated by the 
residual stresses. 

In addition to residual stresses, another factor found to important in accurately predicting 
the initial modulus was proper modeling of the tube cross section. The tubes used in this in- 
vestigation were found to contain resin-rich regions at the ply interfaces. These regions ac- 
counted for seven percent of the overall wall thickness of the tube. The predictions of the 
shear modulus was reduced by ten percent when incorporating the resin "layers" into the 
analysis as well as correcting the material properties for the difference in the fiber volume 
fraction found between the laminated tubes and the panel. 

In order to more accurately predict the shear response, it was found that the tube must 
be modeled using the resin "layers". The inclusion of residual stresses was found to give 
more accurate predictions of tube response in some cases while, in other cases, it did not. 
An offset of the free edge at a longitudinal crack indicated that residual stresses Were present 
in the tubes tested. An improvement can be made to the analytical prediction if an accurate 
determination of the stress free temperature can be made for use In calculating the residual 
stresses. 

The magnitude of the applied stress at failure was found to be path dependent. This de- 
pendence was seen in the tests performed using pure torsion cyclic loading where the tubes 
were subjected to increasingly higher magnitudes of loads in positive and negative torsion. 
Under this sequence of loading, failure occurred at loads lower than had been applied during 
the preceding cycle once a certain stress level was reached. The applied stress at failure, 
then, is a function of the previous loading history. It appears that a critical damage envelope 
(CDE) may exist beyond which failure becomes a path dependent phenomenon. 

Some recommendations for future work in nonlinear analysis of composite tubes are: 

• Incorporate a damage model into the analytical model to predict accumulation of damage 

and its effect on the response of tubes. 
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• Expand the experimental program to Investigate different ratios of torsion and axial loads. 

• Characterize the unloading response of P75/934 graphite/epoxy to verify the assumptions 
made concerning the unloading behavior. 

• Develop model for predicting the magnitudes of the residual stress from the offset of the 
free edge at a longitudinal crack. 

• Perform an experimental investigation similar to that carried out in this research using a 
more nonlinear composite system. 
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Appendix A. Elements of Local Stiffness Matrices 


The radial displacements in the tube were determined by solving a system of equations 
that were generated by assembling a set of local stiffness matrices and force vectors that are 
determined for each layer by writing an expression for the radial component of the radial 
stress at the inside and outside surface of each layer. These expressions are, in matrix form, 
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where a quantity with a k superscript or subscript represents that quantity for the k th layer, a 
superscript represents a quantity at the inside radius of the layer while a "+" superscript 
represents the quantity at the outside radius of the layer. The expressions for the individual 
terms for a transversely isotropic layer are as follows: 
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For layers that are monoclinic or orthotropic, the expressions for the individual terms of 
the local stiffness matrix and force vectors become 
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The integral force equation is of the form 
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and for monoclinic layers, 
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fac? + H?,c? + H?,c? + C? 3 c?)^ — ~ r - k ' ) - (H?,d? + H?,d? + H?,d? + C? 3 d?)f- 
(H|i© k + H?,e? + H?,e? + C? 3 e?)^ ^ - (H?|f? + H?,f? + H^.fg + faf?)^ 
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AT 


F T — 2it f( C 11 a x+ C 12 a 0 + C 13 a r + C 16 a x0)^ * 2 "' ' 

(C^ + ^ 3 )(rJ + ^-W 


L. Ji 

k=1 


c k t 

F, e. 


;k . ~k 


( C 12 + C 13) 


'33 + C 22 

(2fs - - <k + A|,r k-i)( r i< _ Ak - .ta*) 


(I+AkK^-i^) 


,2A V 


(1-W-Cl) 


The integral torque equation is of the form 



k=0 


' F ^N+1 £ o + V{/ N+2Vo~ T x+ F Nl+ F T 


where, for transversely isotropic layers, the coefficients are 


( C 26 ~ £ggKlj r 0 ( C 26 + c 36)( r 1 ~ r o) r Q 

r i + r o 3(r? - rl) 


*k- 

( C 26 + C 3s)( r k ~ r k-l) r k 

( C 26 — c 36) r k r k-1 . 

3(rJ! - rj^,) ( r k + r k-l) 

( C 26 1 ~ C 36 1 ) r k+1 r k ( C 26 1 + C 36 1 )( r k+1 “ r k) r k 

<k+1 + r k 

where k=1,2 ,N-1 

3(<k+1 - 4) 


(C 26 4- C^XrJi — r^,_ 1 )r N 

( C 26 _ < '3«) r N r N-1 

3(r N — r N-i) 

r N + fN-1 
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N 



k=1 


> J'N +2 = 2 ^|( 4 ( iL -^ 1 ) 

k»1 


+ 


/•sk npk 

l .26 ~ 


3Cg2 


(C^ 6 + C3 6 )(r3 — (C; 6 — C^ 6 ) » 3 3 , 

3(4-r 2 k _i) r * + r i<-i 
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‘NL “ “ 27r ^^ j^ P ^ C 36^ 3 k ^ + Y ( C 26 + C 36)[^ r k( ln r k ~ "j" ) ” r k-l( ,n r k-1 “ ’3' ) j * 

k=l 

( C 26 + C 3e)( r k ” r k-l) / 2 . 2 - \ ( C 26 “ ^3e) r k r k-1 , / r k -1 \ 

(rk ln rk " ^ ln ^ H ■ -FT ) 


pk 

‘2 


(C 26 + 2C 36 ) 


(C 2 k 6 + 3C§ 6 ) 


(C 2 e - C 3e) ,24 _ 4 2 

r k + r k _, ' r k r k-i r k r k 


(i^)_ 

;-!>] + P 4 


(C 26 + C 3 e)(r k ~ r k-l) 2 (^26 “ ^3e) ,23 

v k r k-1 


3(^-r 2 k _i) 

\k \,_3 


r k+ r k-1 


r k r k— 1 ) 


(Ci + qW^-rt-^rS-O 


3( r k “ r k-i) 

(Ca e + 4 C^ 6 )^ rk ~ rk ~ 1 ) 


(C 2 e 4- C36)( r k r k -i)( r k ~~ r k -l) (^26 ~ C 36 ) 2 5 


3 (»-k-r k _t) 


(C 26 + 5C 3 k 6 ) 




( C 26 ^36) , 2 6 6 2 

r k + r k _ t ' r k r k-1 ~ r k r k- 


(C 26 + C 36 )(r k f k _i)(r k — r k _^) (C 2e — C^g) 2 7 


r k + rk1 { r k r k-i r k r k— -i) J + 

(C^G + ^k-rtiKrg-rt,) _ 

3( r k — r k _i) 

,) +PS (^6 + 6^6) 


3 (rg - r£_,) 




(H« a ? + ^ 26 a 2 + hW + ^ 3 e a 3 >^ ~ 3 r — ' ■ ) - (Hfebf + H^bJ + + $<$)( ) 


( H ie c i + H 26 c 2 + HgeCg + C 36 c 3 ) 


(Hl 6 e 1 + ^ 26 e 2 + ^ 66 e 6 + C 36 e 3 ) 


(¥*) 


(H k e d k + H 26 d 2 + Hggdg + C 3 gd 3 ) 


(H k e f^ + H 26 f2 + Hggfg + C 36 ^ 


(^) 
* 1(^)1 
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Fj = 2.Y io^16 a x + ^26 a d + ^36 a r + ^66 a x^)f 


-r 3 

r k-1 


AT 


k=1 


\ 


c k T 

F| 6| 


2q 2 

;k , ~k 


/ 3 3 

- k ( r k - r k „ 1 

u 36\ 




3 

3 _3 


( r 2 , n r 4 |„ r ) _ mf ifcl 

3(r?-r? .1 U k k1 k V r k + r k-l V r * 


3(fk — «k-l) 


For monoclinic layers, the terms of the Integral torque equation become 
_ (C20 - A, cleXif Al - rg~ VM' (Css + ^i^3e)( r i + A< - rg + M 

° (2 - — if*0 (2 + ;.,)(rf' - r^') 


'i '„ = 


(C 2 k e + A k C^ + ' , ‘ 


■klM" 


(C 2 6 — ' i k C 3e)( r k ^ 


fk— 1* l ‘) r k lt rk— 1 


(2 + A k )(^-^,) 


(2-2 k )(r^-r 2 +) 


+ 


(C& 1 - - if* (4 t1 + %iCi8 1 )C' 


k+1 


• 2 + ^+i) r jJ k+1 


(2-A k+ ,)(r^r-'k' ,k+ ’) 


.2 A 


(2 + A k+ i)( r kS 


k+1 


- r k 

<f k+ ’) 


where k=1,2,...,N-1 


V F W = 


(C 2 n 6 + A n C 3 n 6 )(^ n + a " - iftftrjr 


( C 26 


2n < '36)0'n ’ l> 




(2 + 2 N )(r 2 -*» 
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/ 3 3 

- k f r k -r kM 


7 =rk 

U 33 ~~ U 22 


f- (C26 + C^6) 


3 3 

'■k”»V_ 1 


(c; e + jj&Xrj* 2 - - rgrtrT * - ■£,*) 
(2 + a k )(r^-r^) 

(5ae ~ ^kCgsjtrg^-r^/- - V k l»,)(^- 2 » - 

(2-2 k )(r^-r 2 k i*i) 


1 + ^k\/i. 2 + Al 


4 4 

'k-'k-'l 


Cg6 — 2 Cqq — k — k / r k — r k _« 

+ #~# (^ 6 + 2 c 5 e)l 

^33 “ u 22 \ 


(Cae + 2 k C^(r 2 k + - i&V (C* 8 - 2 » - r 2+ V^)^- 2 » - ■flfo 

(2 + 2 k )(r 22 « - r|iv (2 - 2 k )(if* - rjft) 


Appendix A. Elements of Local Stiffness Matrices 


162 



Fnl=- 2*) icf, 

Z_j I 


(Cj>e + C3 6 ) 




(Cm + 4<&)(rk* 4 ~ r^, 4 )^ 4 ~ Ci x> ) 
<2 + 4)<r 2 k J "-rft) 


-I f X^ 


(Cm ~ ^C^Ki-j^ft ~ ik + N&Xj 4 ~ fcfr) 
<2-A k )(r^-r&) 


+ G? 


(Cm + 2Cm)( 


(^)- 

(CM + 4^e)(rk +4 -rk^V (Cm - ~ MiWk" *■ - fcfr) 1 


(2 + 2 k )(r^-rft) 


(2-^r^-r^) 


G 3 


(Cm + 3C36> 


(^) 


(Cm + 2 k C$ e )(r£ + 4 - r£, 4 )(r^ 4 -rgl, 4 ) 


(Cm - 2 k CS 6 )(rj;XIi'"‘ - C XlWr *« - ■g? 1 ') 


/-ik \/-2/lt_34- X\m _3f X^JZX^ \/Jl — >t| 


(2 + 2 k )(r^‘-r^) 


.2— /II 


P-rf'-ft 


pk 

g 4 


<cm+«4)( - - / ■ — - ) - 


(Cm + 2 k C§ 6 )(rg ^ - rltHV'-Ch 


(Cm - '’kCMK’-k^Cr 1 - - rr X^)(f k - ik-i'*) 


«-\k \/_24,4+ .4+ At.-2A\, \/_2— 2 


(2 + 4)(r^-r&) 


2 - 


<2-4Xr^-r&) 


pk 

G 5 


/ r k — r k _i ^ 

(C^e - 2 k C38)(r k 4| ' k -i l|t - r k + 4 r k 4 ,)(r k ' 4 - tiN 


(C20 + 6%) I 


(2 + 2 k )(ik 4 — r k 4 i) 


(2-2k)(^‘-r^) 


(Cm + 6Cm) 




(Cm + V&Xrf - rgfrXrg - r£ft) 

(2 + 2k)(r^-r^) 


(Cm - ~ r k +4 r k 4 ,)(rg 4 - 'tt k ) 


(2-A k )<r^-r&) 


(H^a^ + Hm 32 + ^66 a e + cjje a 3) 


(H^c^ + Hijgclj + HggCg + C38C3) 



3 

5 

,5 

«V 

“ r k-1 


5 

' J 

7 

'k 

~ r k-1 


(Hf e b? + Hjgbg + H^bg + C^b^j 


{*£*■)- 

(H^ef + H^eiJ + Hg 8 e£ + ci^e*) ^ k y - - ^ - (H^ff 4- + H^fg + C^fJ) ^ Q k1 ^ j* 


(H^6^ + H^gdjj + Hggdg + ci^d^) 
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( C 16 a x + ^26 a 0 + ^36 a r + C 66 a x6»)f ~ — o'"" 1. 


pk pk ( C 26 + C 36) 
^33 ~ '-'22 


k . pk / rg - \ (C 2 k 6 + A k C k 6 )(rJ + * - r^^)(rg + ^ - fcfy 

26 M 3 ) (2 + ^-r?i V 


(Cm - jiAXjkV^: - r; + S&)(£- ^ - r^Z/*) 


(2-A k )(^-r=i‘ 1 ) 
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Appendix B. Least Squares Polynomial 
Approximation 


For a given layer, the nonlinear strains, eJ* L (r) were approximated by a fifth order 
polynomial to allow for a solution to the governing differential equation for w(r). For a given 
layer and component of strain, the nonlinear strain is approximated by 

e NL (r)-a + br + cr 2 + dr 3 + er 4 + fr 5 

The difference, D, between the true value and the approximation is 

D = a + br + cr 2 *Fdr 3 + er 4 + fr 5 - e NL (r) 

For the least squares method, the squares of this difference is summed over ail of the data 
points within the layer. 


9 - 



i=i 
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where m is the total number of data points. To minimize the function g, the derivative of g 
with respect to each coefficient is set to zero. 


3g ^ 3g 3g 3g 3g 3g 

3a db ~ 3c ~~ 36 ~ de ~~ 3f 


The results of these operations is 


m 

Y^ 


_dg 

da ^ 
i=1 


= 2 ) (a + brj + crf -f- drf -h erf -h frf — 6 j NL ) = 0 


dg 

db 


m 

vn 


= 2 > (a + b^-f erf + drf-f er^ + frf ~^ L )r] = 0 


j=i 


dg 

3c 


m 


2 y (a + brj + erf + drf -f erf + frf - ej NL )rf = 0 


J=i 


dg ^ 
3d 


= 2 y (a + brj + erf + drf + erf + frf - ej NL )r? = 0 


j=i 


m 

- 2^(a + br j + erf + drf + erf + frf - e^ L )rf = 0 
i = 1 
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m 

= 2^(a + br j + cr f + dr f + er t + fr ? _ e j NL ) r f = 0 

j=1 

This system of equations can be solved for the coefficients a, b, c, d, e, and f. in matrix 
form, the system of equations is 


- Zn Z-f Zf Zf Z-f 

a 


V nl 

Z/J 

Z> Z-f Zf Zf Z-f Zf 

b 


2j\ r i 

Z-f Z-f Zf Zf Z-f Z-f 

c 


V NL 2 

2/j r J 

Z-f Z-f Z-f Z-f Z-f Z-f 

d 


V nlj* 
Zj s i r i 

Z-f Z-f Z-f Z-f Z-f Zr 

e 


V N M 
2/J r i 

Z-f Zf Z-f Z-f Z-f Zf 

f 


V NL r 5 

2/j r j 


where each of the summations is from 1 to m. 
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Appendix C. Individual Test Results 
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Figure 48. 0® Tension Response 
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£ x (%) 


Figure 50. 10° Tension Response 
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Figure 52. 10° Shear Response in Tension 
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Figure 53. 45® Tension Response 
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Figure 54. 45* Poisson's Response in Tension 
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Figure 55. 45° Shear Response in Tension 
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Figure 56. 90* Tension Response 


Appendix C. Individual Test Results 


177 




Figure 57. 9Q e Poisson's Response in Tension 
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Figure 58. 0° Compression Response 
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Figure 60. 45° Compression Response 
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Figur* 61. 45* Poisson's Response in Compression 
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Figure 64. 90° Poisson's Response in Compression 
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Figure 65. 0° losipescu Shear Response 
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Figure 66. 90* losipescu Shear Response 
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